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ABSTRACT 


The  Navy  of  tomorrow  will  require  a  robust  and 
reconf igurable  power  system  capable  of  supplying  power  not 
only  to  large  high  power  propulsion  loads  but  to  growing 
combat  system  loads  like  high  power  radar  and  pulse  loads 
such  as  rail  guns  and  f  ree-electron  lasers.  A  critical 
component  in  such  a  system  is  the  phase  controlled 
rectifier.  As  such,  the  issues  associated  with  the 
inclusion  of  a  power  electronics  rectifier  need  to  be 
addressed.  These  issues  include  input  Alternating  Current 
(AC)  interface  requirements,  the  output  Direct  Current  (DC) 
load  profile,  and  overall  stability  in  the  presence  of  non¬ 
linear  loads.  Understanding  these  issues  and  determining 
the  means  of  assuring  compatibility  with  a  Navy  all¬ 
electric  ship  is  the  focus  of  this  thesis. 

By  using  a  Simulink®  model  of  a  variable  parameter 
load,  several  multiple-pulse  count,  high  power  rectifiers 
were  exercised.  The  Simulink®  results  were  compared  to  the 
linearized  small  signal  transfer  function  analysis  results. 

These  experiments  led  to  the  conclusion  that 
increasing  the  pulse  count  and  output  filtering  reduces  the 
input  interface  current  distortion.  However,  there  are 
tradeoffs  in  terms  of  complexity  and  size  of  the  passive 
components,  and  optimization  based  on  source  and  load 
specifications  is  required. 
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EXECUTIVE  SUMMARY 


The  Navy  of  tomorrow  will  require  a  robust  and 
reconf igurable  power  system  capable  of  supplying  power  not 
only  to  large  high  power  propulsion  loads  but  to  growing 
combat  system  loads  like  high  power  radar  and  pulse  loads 
such  as  rail  guns  and  f  ree-electron  lasers.  A  critical 
component  in  such  a  system  is  the  phase  controlled 
rectifier.  As  such,  the  issues  associated  with  the 
inclusion  of  a  power  electronics  rectifier  need  to  be 
addressed.  These  issues  include  input  Alternating  Current 
(AC)  interface  requirements,  the  output  Direct  Current  (DC) 
load  profile,  and  overall  stability  in  the  presence  of  non¬ 
linear  loads.  Understanding  these  issues  and  determining 
the  means  of  assuring  compatibility  with  a  Navy  all¬ 
electric  ship  is  the  focus  of  this  thesis. 

The  objective  of  this  thesis  was  to  investigate 
differences  between  multi-pulse  phase  controlled  rectifiers 
in  continuous  conduction  mode  (CCM)  with  a  two-pole 
inductor-capacitor  (LC)  output  filter.  A  6-pulse  phase 
controlled  rectifier  with  a  two-pole  LC  output  filter  is 
shown  in  Figure  1.  This  circuit  is  the  basic  building  block 
module  for  higher  pulse  count  systems  as  it  is  easily 
stackable . 
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Figure  1 .  Circuit  diagram  of  a  controlled 

three-phase  rectifier  with  a  two-pole 
LC  output  filter.  From  [1] . 


The  key  differences  that  were  evaluated  between  these 
multi-pulse  controlled  rectifiers  were  the  critical 
inductance  Lcrit  (the  inductance  which  is  required  to  maintain 
the  rectifier  in  CCM) ,  the  peak-to-peak  ripple  voltage  as  a 
percentage  of  the  average  output  voltage  %Frc,  /(jx  (also 

commonly  known  as  the  DC  output  voltage  ripple) ,  and  the 
total  harmonic  distortion  in  phase  A  of  the  line  current 
%THDi  . 

The  determination  of  Lcrit  was  made  by  applying  standard 
circuit  analysis  techniques,  first  to  the  6-pulse 
controlled  rectifier  circuit  of  Figure  1,  and  then  by 
extending  the  analysis  to  a  p-pulse  controlled  rectifier.  A 
plot  of  Lcrit  versus  the  firing  angle  a  is  shown  for  a  6-, 
12-,  18-,  24-,  and  96-pulse  controlled  rectifier  in  Figure 
2.  From  Figure  2,  it  is  clear  that  Lcrit  is  lower  for  higher 
pulse  count  controlled  rectifiers  than  for  lower  pulse 
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count  controlled  rectifiers;  the  reduction  in  Lcrit  is  reduced 
by  successively  increasing  pulse  count. 
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Figure  2.  Critical  inductance  plot  for  controlled 

p-pulse  rectifiers  with  unity  resistance 
operating  at  60  Hz. 

The  determination  of  °/oVC(ppiavg)  was  also  made  by  applying 

standard  circuit  analysis  techniques  to  the  6-pulse 
controlled  rectifier  circuit  of  Figure  1  and  extrapolating 
the  results  to  a  p-pulse  controlled  rectifier.  However, 
slight  approximations  were  introduced  in  the  impedance  of 
the  load  and  the  determination  of  the  capacitor  peak-to- 
peak  ripple  voltage.  A  plot  of  0//°^c(pp/ayg)  versus  the  firing 
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angle  a  is  shown  for  a  6-,  12-,  18-,  24-,  and  96-pulse 
controlled  rectifier  in  Figure  3.  From  Figure  3,  it  is 
clear  that  0//°K:(pp/avg)  is  lower  for  higher  pulse  count 
controlled  rectifiers  than  for  lower  pulse  count  controlled 
rectifiers;  the  reduction  in  0//°^c(pp/avg)  is  dramatically 
reduced  by  successively  increasing  pulse  count. 


Figure  3.  Capacitor  peak-to-peak  ripple  voltage  as  a 

percentage  of  the  average  output  voltage  for 
a  6-,  12-,  18-,  24-,  and  96-pulse  controlled 
rectifier  with  a  two-pole  LC  output  filter. 


The  calculation  of  %THDi  was  performed  via  simulation 
using  Simulink®  on  four  p-pulse  controlled  rectifiers  with 
the  same  predefined  passive  components.  A  plot  of  %THDj 
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versus  the  firing  angle  a  for  the  simulated  6-,  12-,  18-, 
and  24-pulse  controlled  rectifiers  is  shown  in  Figure  4. 
From  Figure  4,  it  is  clear  that  %THDi  is  lower  for  higher 
pulse  count  controlled  rectifiers  than  for  lower  pulse 
count  controlled  rectifiers;  the  reduction  in  %THDj  is 
greatly  reduced  by  successively  increasing  pulse  count. 


versus  a  for  p-pulse  controlled  rectifiers 
with  a  two-pole  LC  output  filter. 
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All  of  these  results  led  to  the  main  conclusion  that 
increasing  the  pulse  count  is  beneficial  in  that  it  reduces 
Lcrit ,  °/oVC(pplavgy  and  %THDj  .  However,  in  reality,  tradeoffs 

must  be  considered  in  terms  of  pulse  count  versus 
complexity  and  size  of  the  power  electronic  and  passive 
components . 
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I. 


INTRODUCTION 


A.  BACKGROUND 

A  multi-pulse  phase  controlled  rectifier  is  a  type  of 
alternating  current  (AC)  to  direct  current  (DC)  converter. 
This  means  that  it  takes  an  AC  voltage  source  (typically  a 
sinusoid)  and  changes  it  into  a  DC  voltage.  In  reality,  the 
output  will  be  a  higher  frequency  AC  waveform  with  a 
greatly  reduced  peak-to-peak  voltage  that  has  an  average  DC 
value.  The  term  multi-pulse  usually  applies  to  three-phase 
systems  and  has  to  do  with  the  fact  that  there  are  multiple 
arch-shaped  pulses  corresponding  to  the  peak  regions  of  all 
of  the  different  phase  combinations  of  the  source  voltages 
that  are  seen  by  the  rectifying  circuit  elements.  By  using 
transformers  to  create  phase  shifts  in  the  source  voltages 
and  connecting  these  rectifiers  together  in  different 
fashions,  rectifiers  with  very  high  pulse  counts  can  be 
made.  A  controlled  rectifier  means  that  the  DC  output 
voltage  can  be  controlled  or  adjusted  in  value,  whereas  an 
uncontrolled  rectifier  will  produce  a  fixed  DC  output. 

A  recent  and  growing  trend  toward  an  all-electric  ship 
and  the  development  of  integrated  power  systems  has  caused 
the  role  of  these  types  of  rectifiers  to  become  even  more 
important.  Due  to  the  ever-increasing  shipboard  power 
requirements  to  support  such  items  as  rail-guns,  free- 
electron  lasers,  high  power  radars,  and  electromagnetic 
aircraft  launch  systems,  the  Navy  has  a  tremendous  need  for 
rugged  and  reliable  high  power  conversion  modules  in  the 
megawatt  range  [1]  . 


1 


In  addition  to  the  problem  of  how  to  meet  all  of  the 
power  demands  caused  by  these  new  types  of  loads  is  the 
fact  that  the  power  electronics  (the  phase  controlled 
rectifiers)  supplying  these  loads  are  non-linear.  Non¬ 
linear  loads  produce  unwanted  harmonic  distortion  as  a  by¬ 
product  which  can  negatively  affect  all  directly  connected 
electrical  equipment.  The  phase  controlled  rectifiers  must 
be  capable  of  keeping  the  produced  harmonics  to  an 
acceptable  level  [2] . 

B.  OBJECTIVES 

The  differences  between  multi-pulse  phase  controlled 
rectifiers  in  continuous  conduction  mode  (CCM)  with  a  two- 
pole  inductor-capacitor  (LC)  output  filter  are  investigated 
in  this  research.  The  differences  that  are  specifically 
investigated  are  how  the  critical  inductance,  DC  output 
voltage  ripple,  and  total  harmonic  distortion  (THD)  are 
affected  by  having  differing  pulse  counts. 

C.  RELATED /PREVIOUS  WORK 

In  [3]  and  [4],  many  of  the  equations  for  the  critical 
inductance  of  a  6-pulse  controlled  rectifier  with  a  two- 
pole  LC  output  filter  are  developed  and  presented.  However, 
they  do  not  develop  equations  beyond  those  for  a  6-pulse 
converter.  Similarly,  the  discussion  regarding  DC  output 
voltage  ripple  was  limited  to  that  of  a  6-pulse  converter 
in  [3]  and  [4]  .  However,  there  has  been  much  work  and 
literature  presented  on  the  reduction  of  THD  by  increasing 
pulse  count.  In  [5],  a  24-pulse  count  rectifier  topology  is 
achieved  from  a  single  three-phase  source  using  a  novel 
interconnection  of  conventional  single  and  three-phase 
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transformers  to  effectively  produce  four  three-phase 
systems.  In  [6],  an  18-pulse  rectifier  based  on  using  a 
differential  delta  autotransformer  to  achieve  the  necessary 
phase  shifts  is  presented  for  aerospace  applications.  In 
[7],  a  pulse  doubling  technique  for  a  36-pulse  converter  is 
accomplished  via  two  paralleled  18-pulse  converters  using  a 
delta/polygon-connected  transformer.  In  order  to  double  the 
number  of  pulses  from  36  to  72,  a  tapped  interphase 
transformer  with  two  additional  diodes  is  added  at  the 
output  of  the  rectifier.  In  [8],  the  same  authors  present  a 
similar  pulse  doubling  technique  in  a  12-pulse  converter 
which  uses  two  paralleled  6-pulse  converters  connected 
through  a  star-wired  transformer.  A  tapped  interphase 
reactor  is  connected  at  the  output  to  double  the  pulses 
from  12  to  24.  An  interesting  design  is  presented  in  [9]  by 
way  of  a  parallel-connected  controlled  18-pulse  rectifier. 
Under  reduced  loading  conditions,  the  firing  angle 
increases,  which  results  in  higher  THD  and  lower  system 
efficiency.  To  combat  this  effect,  the  authors  propose  a 
dynamic  control  scheme  which  involves  simultaneous  and/or 
selective  control  of  the  firing  angle  and  converter  pulse 
count.  For  pulse  count  control,  an  appropriate  6-pulse 
rectifier  or  pair  of  6-pulse  rectifiers  is  disabled  or 
enabled  depending  on  the  current  loading. 

D.  THESIS  ORGANIZATION 

Background  material  on  single-phase  rectifiers, 
starting  with  uncontrolled  diode  rectifiers  and  concluding 
with  controlled  rectifiers,  is  provided  in  Chapter  II. 
Critical  inductance  calculations  are  developed,  first  for  a 
6-pulse  controlled  rectifier,  then  for  a  p-pulse  controlled 
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rectifier,  in  Chapter  III.  The  expressions  and  equations 
for  the  DC  output  voltage  ripple,  are  developed  in  Chapter 
IV,  first  for  a  p-pulse  controlled  rectifier  with  a  purely 
resistive  load,  then  for  a  p-pulse  controlled  rectifier 
with  a  two-pole  LC  output  filter.  The  general  concepts  of 
harmonics  and  THD  are  presented  in  Chapter  V,  as  well  as 
how  they  can  be  reduced  by  increasing  pulse  count. 
Simulink®  simulations  of  a  6-,  12-,  18-,  and  24-pulse 
controlled  rectifier  with  a  two-pole  LC  output  filter  are 
provided  in  Chapter  VI.  Conclusions  and  recommendations  for 
future  work  are  discussed  in  Chapter  VII. 
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II.  SINGLE-PHASE  RECTIFIERS 


A.  UNCONTROLLED  SINGLE-PHASE  RECTIFIERS 

1.  Uncontrolled  Single-Phase  Half-Wave  Rectifier 

Before  being  able  to  fully  discuss  what  a  multi-pulse 
phase  controlled  rectifier  in  CCM  with  a  two-pole  LC  output 
filter  is,  it  is  useful  to  begin  with  a  far  simpler 
rectifier  circuit  commonly  known  as  a  half-wave  rectifier. 
One  such  half-wave  rectifier  circuit  is  illustrated  in 
Figure  1,  where  a  sinusoidal  voltage  source  is  applied  to  a 
resistive  load  through  an  ideal  diode.  The  diode  behaves 
like  a  check-valve,  only  allowing  positive  AC  current  to 
pass  through  to  the  resistor.  During  the  negative  portion 
of  the  source  voltage  waveform,  the  diode  blocks  the 
current  flow,  resulting  in  no  voltage  drop  across  the  load. 
Profiles  of  the  source  voltage,  load  voltage,  and  current 
are  shown  in  Figure  2.  From  Figure  2,  it  is  clear  how  the 
circuit  of  Figure  1  earned  its  common  name,  due  to  half  of 
the  source  voltage  input  waveform  being  applied  to  the 
resistive  load  [3],  [4]. 
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2. 


Uncontrolled  Single-Phase  Full-Wave  Rectifier 


A  full-wave  rectifier  is  the  next  logical  progression 
beyond  a  half-wave  rectifier.  A  simple  full-wave  rectifier 
(sometimes  called  a  bridge  rectifier)  with  a  resistive  load 
is  shown  in  Figure  3. 


<L 


Figure  3.  Uncontrolled  full-wave  rectifier  circuit 

with  ideal  diodes  and  a  resistive  load.  From  [3] . 

In  this  circuit,  the  source  voltage  and  load  resistor 
are  identical  to  those  of  the  half-wave  rectifier  case.  In 
the  full-wave  rectifier  circuit,  in  addition  to  the  current 
that  flows  during  the  positive  portion  of  the  source 

voltage  waveform,  current  will  also  flow  during  the 

negative  portion  due  to  the  second  pair  of  diodes.  During 
the  positive  portion  of  the  source  voltage  waveform,  diodes 
D,  and  D4  are  "on"  while  diodes  D1  and  Ds  are  "off, "  and 
current  will  flow  in  the  direction  indicated  by  the  load 

current  iL  .  During  the  negative  portion  of  the  source 

voltage  waveform,  diodes  D{  and  D4  are  "off"  while  diodes 
D1  and  D 3  are  "on."  Notice  that  the  direction  of  iL  is 
unchanged.  This  unidirectional  current  results  in  a 
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positive  load  voltage  drop  during  both  the  positive  and 
negative  portions  of  the  source  voltage  waveform.  Profiles 
of  the  source  voltage,  source  current,  load  voltage,  and 
load  current  as  described  are  illustrated  in  Figure  4  [3], 

[4]  . 


Figure  4.  Source  voltage,  source  current,  load 

voltage,  and  load  current  waveform  profiles  for  an 
uncontrolled  full-wave  rectifier  circuit.  After  [3]. 


B.  CONTROLLED  SINGLE-PHASE  RECTIFIERS 

1.  Controlled  Single-Phase  Half-Wave  Rectifier 

If  the  diode  from  the  circuit  of  Figure  1  is  replaced 
with  a  thyristor  or  silicon  controlled  rectifier  (SCR) ,  the 
result  is  a  controlled  half-wave  rectifier  of  Figure  5.  The 
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SCR  behaves  similarly  to  a  diode  in  that  it  is  a  one-way 
device  and  will  block  current  flow  in  the  negative 
direction.  However,  it  will  not  conduct  in  the  forward 
direction  until  an  appropriate  trigger  signal  has  been 
applied  to  its  gate  [10] .  In  reality,  there  will  be  a  short 
delay  before  the  SCR  turns  "on"  even  after  it  has  been 
adequately  triggered;  but  for  the  purpose  of  explanation, 
the  SCR  in  Figure  5  can  be  considered  ideal  and  will  turn 
"on"  immediately  when  triggered.  The  angle  at  which  the  SCR 
is  triggered  is  commonly  called  the  firing  angle  a  .  If  the 
SCR  in  Figure  5  is  triggered  at  time  tx  ,  which  is  equal  to 

a  I  co  (where  co  is  the  radian  frequency),  the  resistor 
voltage  and  current  will  be  as  shown  in  Figure  6  [3],  [4]  . 


Figure  5.  Controlled  half-wave  rectifier  circuit 

with  an  ideal  SCR  and  a  resistive  load.  From  [3] . 
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1 1  T/2  T  t 

Figure  6.  Source  voltage,  load  voltage,  and  load 

current  waveform  profiles  for  a  controlled  half-wave 
rectifier  circuit  triggered  at  time  tx  .  From  [3]  . 

2.  Controlled  Single-Phase  Full-Wave  Rectifier 

By  replacing  the  diodes  in  Figure  3  with  SCRs,  a 
controlled  full-wave  rectifier  is  created  as  shown  in 
Figure  7.  SCRs  are  gated  in  pairs  with  the  firing  angle  a 
again  corresponding  to  time  tx=alco.  During  the  positive 
portion  of  the  source  voltage  waveform,  all  SCRs  are 
initially  "off"  until  SCRi  and  SCR4  are  gated  "on"  while 
SCR2  and  SCR3  remain  "off"  due  to  being  reverse  biased. 
During  the  negative  portion  of  the  source  voltage  waveform, 
the  opposite  situation  occurs  and  all  SCRs  are  "off"  until 
SCR2  and  SCR3  are  gated  "on"  while  SCRi  and  SCR4  remain 
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reverse  biased  and  "off."  The  resulting  current  yields  zero 
voltage  drop  across  the  load  resistor  during  the  time  the 
SCRs  are  "off"  and  a  positive  load  voltage  drop  during  both 
the  positive  and  negative  portions  of  the  source  voltage 
waveform  when  the  appropriate  SCRs  have  been  gated  "on." 
Profiles  of  the  source  voltage,  source  current,  load 
voltage,  and  load  current  as  just  described  are  illustrated 
in  Figure  8  [ 3 ] ,  [  4 ]  . 
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Figure  7.  Controlled  full-wave  rectifier  circuit 

with  ideal  SCRs  and  a  resistive  load.  From  [3] . 
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Figure  8.  Source  voltage,  source  current,  load 

voltage,  and  load  current  waveform  profiles  for  a 
controlled  full-wave  rectifier  circuit  triggered  at 
time  tx  .  t2=tx+T  /  2.  After  [3]. 
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III.  CRITICAL  INDUCTANCE 


A  controlled  three-phase  rectifier  with  a  two-pole  LC 
output  filter  is  depicted  in  the  bridge  configuration  of 
Figure  9.  In  this  configuration,  the  number  of  pulses  is 
twice  the  number  of  phases,  hence,  three-phase  becomes 
synonymous  with  6-pulse  [11]  .  One  very  useful  circuit 
element  parameter  for  these  types  of  rectifiers  is  the 
value  of  the  inductor  necessary  to  maintain  the  circuit 
operating  in  COM.  This  particular  value  of  inductance  is 
called  the  critical  inductance  Lcrit  .  For  nominal  circuit 
operation,  the  instantaneous  inductor  current  is  required 
to  always  be  positive,  which  is  the  definition  of  COM,  and 
the  capacitor  is  required  to  be  very  large  such  that  the 
instantaneous  capacitor  voltage  is  approximately  equal  to 
the  average  capacitor  voltage  vc(t)  =  FC(avg)  [3],  [4].  While  it 

does  require  lengthy  equations  and  calculations,  we  can 
solve  for  Lcrit ,  as  will  be  shown  in  the  continuing  sections. 
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filter .  From  [ 4 ] . 
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A. 


AVERAGE  LOAD  VOLTAGE  CALCULATION 


To  aid  in  understanding  the  relationship  between  the 
voltage  sources  for  a  three-phase  circuit  such  as  the  one 
shown  in  Figure  9,  the  three  line-to-line  source  voltages 
are  plotted  in  Figure  10.  While  each  voltage  waveform 
possesses  the  same  magnitude,  they  are  offset  from  each 
other  by  120 ° . 


Figure  10.  Line-to-line  source  voltage  waveforms 

Vah  r  Vbc ,  and  Vca  with  unity  magnitudes  for  a 
hypothetical  three-phase  circuit. 

In  addition  to  the  three  line-to-line  source  voltages 

as  indicated  by  the  polarities  in  Figure  9,  there  are  three 

more  line-to-line  voltages  that  exist  by  looking  at  them 

14 


Voltage  (Volts) 


from  the  opposite  polarities.  These  additional  voltage 
waveforms,  along  with  the  originals,  are  plotted  in  Figure 
11 . 


Figure 


11 . 
v 

v  cb  ' 


magnitudes 


Line-to-line  voltage  waveforms 

Vab’  VaO  Vbc>  Vba'  and  Ka  with  Unity 

for  a  hypothetical  three-phase  circuit. 


From  Figure  9  and  [3]  and  [4],  the  average  capacitor 
voltage  Vc(avg)  is  given  by 


a+ 120° 


Va°vg)  =  7^  I  Kn sin  ydr  < 


a+ 60° 


(l) 
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where  Vm  is  the  peak  value  of  the  line-to-line  voltage.  The 
six  in  the  numerator  is  due  to  the  fact  that  there  are  six 
voltage  pulses  over  one  period  for  a  three-phase  rectifier, 
as  illustrated  in  Figure  11.  The  integration  limits  are 
solved  for  by  first  dividing  one  period  of  360°  by  six 
pulses  to  obtain  a  pulse-width  of  60°.  Dividing  this  pulse- 
width  by  two  and  adding  and  subtracting  it  with  90°  yields 
the  upper  and  lower  limit  firing  angle  offsets  for  the 
upper  and  lower  integration  limits.  By  integrating  (1),  we 
get 

v 

’  C(avg) 

Evaluating  (2),  we  obtain 

W 

VC(avg)  = - — [cos(or  + 120°)  -  cos(a  +  60°)]  ,  ( 3 ) 

n 

which,  by  applying  the  trigonometric  identity, 
cos(a  +  b)  -  cos(fl)cos(6)-sin(fl)sin(&) ,  becomes 


3F 


n 


-cos/ 


v=rY+hQ° 


(2: 


3(/  fcos(a)cos(1200)-sin(a)sin(120°)  j 
(avs)  n  |-[cos(a)cos(600)-sin(a)sin(60°)]j  ' 


which  ultimately  simplifies  to 

3V, 


V, 


C(avg) 


K 


-cos  a 


(4) 


(5) 


B.  PROCEDURE  TO  SOLVE  FOR  THE  BOUNDARY  ANGLE  ah  AS  THE 
TRANSITION  POINT  BETWEEN  SUB-MODES  OF  CCM 

In  the  circuit  of  Figure  9  there  exists  a  firing  angle 
ah  such  that  the  instantaneous  source  voltage  is  egual  to 

the  average  capacitor  voltage  |v5(a)| -VC{avg)  .  Therefore,  for  an 
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as  v, 


instantaneous  source  voltage  vs(t)  given  as  v^(t)  =  Vms\n(cot)  and 

an  average  capacitor  voltage  previously  solved  for  in  (5)  , 
we  have  the  equality 


3V 

Vm  sin(tftfj)  =  — -cos ab 
n 


(6; 


By  substituting  a>tx  =  ah  +  60°  ,  we  have 


3V„, 


Vm  sin(«A  +  60°)  =  — -cos  a, 


n 


By  applying  the  trigonometric  identity 

sin(a  +  b)  =  sin(<7)cos(&)  +  cos(o)sin(&) ,  we  produce 


(7: 


3V„ 


Vm  [sin ( ah )  cos  ( 60°)  +  cos  ( ah )  sin ( 60°)]  =  — —  cos  ah , 


(8! 


which  can  be  rewritten  as 


lr  . 


2  L 


sin  ah  +  yf3  cos  ah 


=  —  cos«A 
n 


(9; 


Grouping  similar  terms  yields 


sin  ab  = 


f— -  V3 

\n 


cos  a.  , 


which  can  then  be  reduced,  leaving 


tan  ab  = 


(io; 


(ii) 


which  can  be  solved  for  ah  to  two  decimal  places  as  [3],  [4] 


a,  =  10.08°  . 


(12: 
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C.  PROCEDURE  TO  SOLVE  FOR  THE  CRITICAL  INDUCTANCE  FOR  A 

FIRING  ANGLE  >  10.08° 

With  a  continuous  inductor  current  and  for  firing 
angles  greater  than  10.08°,  the  absolute  value  of  the 
instantaneous  source  voltage  is  greater  than  the  average 
capacitor  voltage  at  the  point  of  commutation,  given  by 
«  +  60°  =  &tf1.  Therefore,  the  slope  of  the  inductor  current  at 
the  point  of  commutation  is  always  greater  than  zero,  and 
it  can  be  said  that  the  voltage  source  is  always 
electrically  connected  to  the  load  via  the  SCR  bridge 
network.  In  Figure  12,  the  waveform  profiles  of  the 
instantaneous  load  voltage,  average  load  voltage,  and 
instantaneous  inductor  current  for  the  circuit  of  Figure  9 
are  provided.  In  order  to  establish  a  well-defined  boundary 
condition  for  the  proceeding  calculations,  we  let  L  =  Lcrjt  so 
that  iL(tl)  =  iL(t2)  =  0  every  one-sixth  period  of  the  source 
voltage,  resulting  in  operation  on  the  boundary  between 
continuous  and  discontinuous  modes  (also  known  as  being 
barely  continuous),  as  shown  in  Figure  12. 
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f 


2 


Figure  12.  Load  voltage,  average  load  voltage. 


and  inductor  current  waveform  profiles  for  a 
controlled  three-phase  rectifier  with  a  two-pole  LC 
output  filter  with  a  firing  angle  greater  than  10.08°. 

After  [ 3 ] . 

From  Kirchhoff '  s  Voltage  Law  (KVL) ,  we  have 


(13) 


where  vs(t)  is  the  source  voltage  given  by  vs(t)  =  Vmsm{cot)  . 


Substituting  (5)  for  VC(avg)  and  L  for  vx(t)  ,  we  get 


(14) 


Rearranging  (14),  we  get 


y  • 


(15) 
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To  solve  for  iL(t) ,  it  is  necessary  to  set  up  the  integral  of 
(15) ,  leading  to 


V  r 


sin(fttf) - cos  a 

n 


dt , 


(16; 


which,  when  integrated,  results  in 


V 

L  L 


-cos  (cot)  3tcosa 


co 


n 


Jr,  * 


(17; 


Upon  evaluating  (17)  and  factoring  out  1  /  co  ,  we  get 


V 

iL(t)  = ^7 

coL 


-cos  {cOt)~ 


3cotcosa 


n 


3  cot,  cos  a 


-cos(ft>^) - L 

v  n  J 


(18; 


Rearranging  (18),  we  get 

V, 


4(0 =-7 

COL 


f  -i  5 


cos(  cot  -  cos  (cot)  + j[cot{  -  cot)  cos  a 


(19; 


By  substituting  <^=<2  +  60°  in  (19),  we  get 


V 

zi(0  =  — y 

coL 


cos(a  +  60o)-cos((«0  + 


i  3  V  an  n  ^ 
+ - cot 


\n  j 


180  3 


COSCf 


(20 ; 


It  is  now  necessary  to  calculate  the  average  load  current 
Javg  >  given  by 


I avg  rp 


1  T 

-\iL{t)dt 


(21 : 


Substituting  (20)  into  (21),  we  get 


I  ~  6V” 
avg  IncoL 


a+ 120° 


a+6  0° 


cos(a  +  60°)  -  cos(cot)  +  — 

n 


3  f  an  n  ^ 
■  + - cot 


180  3 


cos  a 


d(cot)  .  (22 ) 
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Integrating  and  evaluating  the  integral  given  by  (22)  in 
one  step  would  be  very  cumbersome,  so  it  is  broken  up  into 
three  parts.  Let 


<*+120° 

Xl  =  cos(a  +  60°)  J  d(cot) , 

a+ 60° 

a+ 120° 

X2  =  J  \-cos(a>t)\l((ot) , 

a+ 60° 


and 


X3  = 


3  cos  a 


K 


a+ 120° 


a+ 60° 


an  n 

- + - cot 

180  3 


d(cot) 


Starting  with  Xx  and  integrating  (23) ,  we  get 


[/  ,^\-ifi#=a+120° 

cot  cosy  a +  60' )  , 

v  >  -lwf=a+60° 

which,  when  evaluated,  becomes 

Xx  =ycos(a  +  60°)  . 

By  applying  the  trigonometric  identity 
cos(a  +  6)  =  cos(tf)cos(6)-sin(a)sin(6)  to  (27),  we  get 


Xj  =y[cos(a)cos(60°)-sin(a)sin(60°)] , 


which  can  be  rewritten  as 


X,  = 


n 


cos(a)- V3sin(a) 


Continuing  with  X2  and  integrating  (24),  we  get 


Tjr  V  -  /  \  ~\cot=a+\20° 

X2  =  [-sin(e#)]  , 


which,  when  evaluated,  becomes 
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(23) 


(24) 


(25) 


(26) 


(27) 


(28) 


(29) 


(30) 


X2  =sin(a  +  60°)-sin(a  +  120°)  . 


(3i : 


By  applying  the  trigonometric  identity 
sin(a  +  b)  =  sin(tf)cos(&)  +  cos(tf)sin(&)  to  (31),  we  get 


X2  =  sin(«)  cos(60°)  +  cos(a)  sin(60°) 
-[sin(<2)cos(120°)  +  cos(«)sin(120°)] 


(32: 


which  simplifies  to 


X2  =  sin  a 


(33) 


Finishing  with  X2  and  integrating  (25) ,  we  get 


^3  = 


3  cos  a 

V  n  J 


f  an  n  ^ 
-  +  — 


180  3 


cot- 


(coty 


-\cot=a+\20o 


-I  cot=a+6  0° 


which,  when  evaluated,  becomes 


(34; 


X3  = 


3  cos  a 
V  n  J 


f  an  n\(  n\  1  I  (  an  2  n 
+  — 


180  3 


n 

v  3  j 


v 


180  +  3  j 


an  n 
■  +  — 


180  3 


os; 


Simplifying  (35) ,  we  get  the  result 


^3=- 


;rcos<2 


(36; 


Inserting  the  resultant  expressions  from  (29),  (33),  and 

(36)  back  into  (22),  we  have 


n 


l  =  6V"‘  j 

avg  2ncoL  1  6 


cos(a)- V3sin(a) 


+  sina-- 


^■cosa 


(37: 


which,  when  simplified,  produces 

3V„,  sin  a 


I avg 


ncoL 


os: 
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in  (5)  ,  we  have 


Because  we  have  already  solved  for  FC(flVg) 
another  expression  for  /  using  Ohm' s  Law,  as 


avg 


— -cosa  . 
7lR 


(39) 


By  substituting  (39)  into  (38),  we  can  solve  for  Lcrit  as  [3], 
[4] 


^ crit 


Rian  a 


co 


n 


■S 


(40) 


D.  PROCEDURE  TO  SOLVE  FOR  THE  CRITICAL  INDUCTANCE  FOR  A 

FIRING  ANGLE  <  10.08° 

With  a  continuous  inductor  current  and  for  firing 
angles  less  than  10.08°,  the  absolute  value  of  the 
instantaneous  source  voltage  is  less  than  the  average 
capacitor  voltage  at  the  point  of  commutation.  Therefore, 
the  slope  of  the  inductor  current  at  the  point  of 
commutation  is  always  less  than  zero.  Unlike  the  previous 
sub-mode  in  which  the  minimum  inductor  current  occurred  at 
the  firing  point,  for  this  sub-mode,  the  minimum  inductor 
current  occurs  at  an  angle  d^  =  coti  (where  6  =  #3-60°  referenced 
to  a  ) ,  as  shown  in  Figure  13. 
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Figure  13.  Load  voltage,  average  load  voltage,  and 
inductor  current  waveform  profiles  for  a 
controlled  three-phase  rectifier  with  a  two-pole  LC 
output  filter  with  a  firing  angle  less  than  10.08°. 

After  [ 3 ] . 


Once  again,  we 

can  say  that 

the 

voltage 

source  is 

always 

electrically 

connected  to 

the 

load  via 

the  SCR 

bridge 

network,  and 

the  conducting 

SCRs 

are  able 

to  ride- 

through  the  minimum  current  point  at  t3  .  At  t3 ,  the  average 
capacitor  voltage  is  equal  to  the  instantaneous  source 
voltage  and  the  slope  of  the  inductor  current  is  zero. 
Therefore,  for  an  instantaneous  source  voltage  vs(7)  given 

as  vs  (t)  =  Vm  sin(&tf3)  and  an  average  capacitor  voltage  previously 
solved  for  in  (5) ,  we  have  the  equality 

3V 

Fmsin  (cot3)  =  — —  cos  a  .  (41) 

n 

By  substituting  <x>t3=d  +  60°,  we  have 
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Vm  sin(#  +  60°) 


(42) 


3V„, 


-cos a  , 


n 


which  can  be  solved  for  9  as 


9  =  sin  1 


—  cos  a 
\n  j 


-60c 


(43) 


For  t>t3,  the  instantaneous  source  voltage  is  greater  than 
the  average  capacitor  voltage,  and  the  inductor  current 
begins  to  increase.  By  triggering  the  SCRs  at  ah  as 
previously  determined  in  (12),  an  additional  verification 
can  be  made  at  the  boundary  between  sub-modes.  For  the 

instantaneous  source  voltage  vs(0  -  Ki  sin(W3)  and  an  average 
capacitor  voltage  as  previously  solved  for  in  (5)  (but  with 
ab  now  substituted  for  a )  ,  we  have  the  equality 


3V 

Vm  sin  ( a>t3 )  =  — -  cos  ah  . 
n 

By  substituting  coti=9h  +  6 0° ,  we  have 


(44) 


3V„, 


Vmsm(9b+60°)  =  — -cos ab , 


n 


which  can  be  solved  for  6b  as 


9h  =  sin  1 


f 


—cos  ah 

\7T 


-60°  . 


(45) 


(46) 


At  ab- 10.08°,  9h  =  10.08°.  In  order  to  once  again  establish  a 
well-defined  boundary  condition  for  the  subsequent 
calculations,  we  let  L  =  Lcrit  so  that  iL(t3)  =  0  every  half-period 
of  each  phase  of  the  source  voltage,  again  resulting  in 
operation  on  the  boundary  between  continuous  and 
discontinuous  modes,  as  illustrated  in  Figure  13.  We  begin 
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once  again  with  the  KVL  equation  as  given  in  (13),  with  the 
same  expressions  for  vs(/),  vt(7)  ,  and  VC(avg) ,  working  through 
identical  equations  given  by  (14)  and  (15),  and  arrive  at 


4(0 


K 

L 


3 

- cos  a 

n 


(47) 


After  integration,  we  have 

V 

4(0  =  — 

L 

and  after  evaluating  (48)  and  factoring  out  H  co ,  we  get  the 
result 


-cos(&>h)  3tcosa 


co 


n 


(48; 


V 

4(0  =  -7 

coL 


-cos(cot)- 


3cotcosa 


K 


-cos  (cot3)~ 


3  cot 3  cos  a 
n  ) 


(49; 


Rearranging  (49),  we  get 

V„ 


4(0  =  -7 

coL 


^  3  ^ 


cos(ft#3 )  -  cosfwh)  +  —  ( cot3  -  cot )  cos  a 

7 


:  50 : 


By  substituting  <yf3=6,  +  60°  into  (50),  we  get 


F 

4(0  =  — 7 

coL 


cos  ( 0  +  60°)  -  cos(d»0  + 


7-2  1 


On  n 


-  + - cot 


180  3 


cos  a 


(51 : 


Once  again,  we  calculate  the  average  load  current  using 
(21).  By  substituting  (51)  into  (21),  we  get 


I 

avg 


6V 

_ m 

2  nco  L 


a+ 120° 


a+6  0° 


cos  (6*  +  60°)  -  cos  (cot)  +  — 

71 


f  On  n 

- + - cot 

v180  3 


cosa 


d(cot) 


(52) 


Integrating  and  evaluating  the  integral  given  by  (52)  in 
one  step  would  once  again  be  very  cumbersome,  so  it  is 
separated  into  three  parts.  Let 
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a+ 120° 

X1  =  cos(<9  +  60°)  J  d(cot) , 

a+6  0° 


(53) 


X2=X2  from  (24),  and 


^3  = 


3  cos  a 


n 


a+ 120° 


a+60° 


-  +  —  cot 


180  3 


d(cot) 


(54: 


Starting  with  Xx  and  integrating  (53) ,  we  get 


Xx  -  [<^cos(0  +  6O0)] 


cot=a+ 120° 
6tf=a+60° 


(55: 


which,  when  evaluated,  becomes 


Xx  =  ycos(6,  +  60°) 


(56; 


Since  X2  is  identical  to  the  previous  X2  ,  the  result  is 
provided  from  (33)  .  Finishing  with  X3  and  integrating  (54), 
we  get 

cot=a+ 120° 

,  (57) 

cot=a+ 60° 


^3  = 


^3cos 

V  71  J 


6k  K 


180  3 


cot- 


(coty 


which,  when  evaluated,  becomes 


^3  = 


3  cos  o' 


V  n  ) 


6k  n 
180  Jj 


54^-3 


1  If  CCK  2k 

vl80  T 


\2  4  \2 

'  '  a,K  k  ' 

v180  y 


1  j 


(58: 


Simplifying  (58),  we  get 


X3  = 


^  n  5 


Vl80y 


(6,-a-30°)cosa  . 


(59) 


Inserting  the  resultant  expressions  from  (56) ,  (33) ,  and 

(59)  back  into  (52),  we  have 
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(60) 


3F  jt  (  jt  \ 

Inv„  = — —  —  cos  ( 0 +  60°)  + sin  a  +  (0  -  a -30°)cosa 

avs  ttcoL  |_  3  V  7  U80jV  7 

Using  (39)  once  again,  we  can  substitute  it  into  (60)  and 
solve  for  Lcrjt  as 

L.= — — —  —  cos(#  +  60°)  +  sina  +  — —  (0-a -30°)cosa  (61) 

wcosa |_ 3  ^180 y  J  . 

Although  61  is  a  function  of  a  as  given  by  (43)  , 
substituting  it  into  (61)  would  only  serve  to  make  it  more 
unwieldy,  so  we  leave  it  in  its  present  form  [3],  [4] . 

E.  SUMMARY  OF  CALCULATIONS  FOR  Lcrit 

To  summarize,  we  have  the  following  equations  for  Lcrit  : 


R  tan  a  n  V3 

®  l  6  , 


for  a>  10.08° 


(62) 


From  Figure  14,  the  critical  inductance  for  the  6-pulse 
rectifier  circuit  of  Figure  9  can  be  seen,  assuming  a  unity 
resistance  and  a  frequency  of  60  Hz.  An  important 
conclusion  is  that  as  a  increases,  Lcrit  increases  at  an 
increasing  rate.  This  is  the  key  reason  that  circuit 
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designers  of  controlled  rectifiers  desire  to  keep  a 
relatively  small— the  size,  weight,  and  cost  of  the  inductor 
cannot  be  ignored. 


x  10'3 


0  10  20  30  40  50  60  70  80  90 

a  (degrees) 

Figure  14.  Critical  inductance  plot  for  a  controlled 
6-pulse  rectifier  with  unity  resistance 
operating  at  60  Hz. 

F.  CRITICAL  INDUCTANCE  FOR  A  CONTROLLED  P-PULSE  RECTIFIER 

One  way  to  increase  the  pulse  count  of  this  type  of 
rectifier  is  by  stacking  6-pulse  rectifiers  in  series— for 
instance,  two  6-pulse  rectifiers  could  be  connected 

together  to  make  one  12-pulse  rectifier.  While  this  has 
several  benefits,  an  important  one  is  that  it  decreases  the 


29 


value  of  Lcrit  .  To  prove  this,  one  could  certainly  go  back  and 
tediously  derive  the  necessary  equations,  solve  for  Lcrit ,  and 
plot  it.  Indeed,  that  is  exactly  what  was  done  for  the  12-, 
18-,  and  24-pulse  rectifier  circuit  configurations,  the 
workings  of  which  have  been  spared  from  this  thesis.  The 
purpose  of  this  labor  was  to  uncover  a  possible  pattern. 


This 

work 

was  not  done 

i  n  va  i  n  ; 

a  clear  pattern 

exists . 

There 

are 

equations  to 

calculate 

Lcrit  for  a  controlled  p- 

pulse 

rectifier.  These 

equations 

are  derived  in 

the  same 

manner  as 

before,  only  this 

time  with  some 

clever 

substitutions.  First,  let  p  equal  the  pulse  count.  Then, 
let  the  lower  integration  limit  0L  be  expressed  in  degrees 
as 


e,  =  90° 


"l80°" 

V  P  V 


(65) 


and  the  upper  integration  limit  0V  ,  also  in  degrees,  be 

180° 


0U  =  90°  + 


V  P  J 


(66; 


thereby,  making  both  functions  of  pulse  count.  Next,  let 
A  =  cos 0L  ,  B  =  sin 0L  ,  C  =  cos0ur  and  D  =  sin0LI  ;  therefore,  C  = -A 
and  D  —  B  . 


1 .  Average  Load  Voltage  Calculation 

With  our  new  substitutions,  the  expression  for  the 
average  load  voltage  from  (1)  now  becomes 

CC+Qy 

Vc(avg)  =  j  vm  sin /dr,  (67) 

lK  A 
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which,  when  integrated,  yields 


vc(aVg)  =-^pL[cos(a  +  0u)-cos(a  +  0L)]  . 
2.71 


(68; 


By  applying  cos(a +  b)  -  cos(a)cos(b)-sin(a)sm(b)  to  (68),  we  get 


V, 


C(avg) 


pV 

-  — -2l[c cos(a)  -Dsin(a)  -Acos(a)  +  f?sin(a)] , 


2  n 


(69; 


which  simplifies  to 


Vc(avg)=^JJL™sa 

n 


(7o; 


2 .  Procedure  to  Solve  for  the  Boundary  Angle  ab  as 
the  Transition  Point  Between  Sub-modes  of  CCM 

By  beginning  with  (6),  but  this  time  substituting 
a>t1=ah  +  0L  and  our  new  expression  for  average  capacitor 
voltage  from  (70),  we  have 

Vm  sin  (ay  +  0L)=  pVm  cos  ab  .  (71) 

n 

By  applying  the  trigonometric  identity 

sin(a  +  b)  =  sin(tf)cos(&)  +  cos(a)sin(Z?) ,  we  produce 

Vm  [sin (ab) cos (0L)  +  cos (a6)sin(0L)]  =  ^^Lcosab ,  (72) 

which  can  be  rewritten  as 


A  n 

^sin(a6)  +  5cos(ai)  =  — -cosor 


Grouping  similar  terms,  yields 


As\n(ab)  = 


r  ±-B) 

V  n  j 


cos  a,  , 


(73) 


(74) 
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which  can  then  be  reduced,  leaving 


,  n.  ( p  B 

tan  (ab)  = - - 

V  71  A. 


Therefore,  we  have 


_i  (  p  B 

a,  =  tan  - - 

b  l  n  A 


3 .  Procedure  to  Solve  for  the  Critical  Inductance 
for  a  Firing  Angle  >  ab 

Starting  with  (13),  we  use  the  same  original 
substitutions,  with  the  exception  of  our  new  expression  for 
VC(avg) ,  and  obtain 

Vm  sin(ft#)  =  L  —  +  CQS a  ' 
dt  n 

which,  after  we  solve  for  diL/dt,  integrate,  and  simplify  as 
before,  we  get 


V  (Ap  | 

iL(t)  =  —SL-  cos(ojtt  )  -  cos(ojt)  +  —  (co^- cot)  cos  a 
coL  \  x 


By  substituting  (ot{=a  +  dL  in  (78),  we  get 


iL(t)  =  — —  cos(a  +  0L)  -  cos(cot)  + 


Ap  V  an  0Ln 


n  A 180  180 


+  — - cot  cosa 


Calculating  I  by  substituting  (79)  into  (21),  we  get 


a+du  cos{a  +  6L)-cos(oot) 

hvg  =  —j  j  ( Ap  V  an  0Ln  \  d{cot) 

2ncoL  *  +  — -  - +  — - cot  cos  a 

L  l  n  180  180 
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Breaking  up  the  integral  into  similar  parts  as  before,  we 
have : 


W.  I  VfJ 

Xx  =  cos(a  +  0L)  |  d(cot) , 


X2  =  J  [~cos(cot)]d(cot) , 


Apcosa  f '[  an  6,n 

X,=— -  - +— - cot  d( cot) 

n  All  80  180  / 


Starting  with  Xx  and  integrating  (81),  we  get 


Xx  =  \jotcos(Ka  +  0L) 


cot=a+Ou 


which,  when  evaluated,  becomes 


X1=  —  (3/-0i)cos (a  +  0L)  . 

V  ioU  / 


We  apply  cos(a  +  b)  =  cos(tf)cos(&)-sin(tf)sin(&)  to  (85),  yielding 


xi  =  —  (^/-^)[cos(a)cos(^)-sin(a)sin(^)]  , 
v  loU  y 


which  can  be  rewritten  as 


Xx  =  —  r^lcos(a)-JSsin(a)l 
l  P  ) 


Continuing  with  X2  and  integrating  (82),  we  get 


v=[-Sin 


which,  when  evaluated,  becomes 
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X2  =sin(a  +  6,i)-sin(a  +  6,(7)  .  (8  9; 

By  applying  sin(a +  b) -sm(a)cos(b)  +  cos[a)sin(b)  to  (89),  we  get 


X2  =  sin(a)cos(6,L)  +  cos(<2)sin(6,i) 
-[sin(a)cos(6,L,)  +  cos(a)sin(6,(/)] 


which  simplifies  to 


X2  =  2  A  sin  a 


Finishing  with  X2  and  integrating  (83),  we  get 


X3  = 


Ap  cos  a 
\  *  ) 


an  G.n 
■  +  ■  L 


180  180 


^  (cot) 

I  cot-- — — 

J  2 


cot=a+0T, 


cot=a+0. 


which,  when  evaluated,  becomes 


X3  = 


Ap  cos  a 

\  n  j 


( an  0,n^ 

\i0u-eA 

(  n  ^ 

- +  — 

OO 

o 

00 

o 

|  V  u  L  '  \ 

U8oJ 

1  I  (an  0„n\  (  an  6,n 

■+  L 


180  + 180 


180  180 


Simplifying  (93),  the  result  is, 

_  2Ancosa 
P 


r  _  PVm 
avs  2ncoL 


/  2n^ 
l  P  j 


[Acosa  -5  sin  a]  +  2^4  sin  a  - 


2  An  cos  a 


which,  when  simplified,  produces 

V„  sin  a 


“vg  ncoL 


( Ap-nB )  . 


(9o; 


(9i: 


(92: 


(93; 


(94; 


Inserting  the  resultant  expressions  from  (87),  (91),  and 

(94)  back  into  (80),  we  have 


(95; 


(96; 
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in  (70)  ,  we  have 


Because  we  have  already  solved  for  Fc(flVg) 
another  expression  for  Iavg  using  Ohm' s  Law,  as 


,  ApVm 
Iavg= — -f-cosa  . 

7TK 


(97: 


By  substituting  (97)  into  (96),  we  can  solve  for  Lcrit  as 

f  — 


L  .  = 

cnt 


Rian  a 


co 


1- 


7lB 

Ap 


os: 


4 .  Procedure  to  Solve  for  the  Critical  Inductance 
for  a  Firing  Angle  <  ab 

Starting  with  (41),  we  substitute  our  new  expressions 
for  VC(avg)  and  cot3  to  obtain 


Vm  sin  (0  +  0L)  =  — ^^-cos  a 

K 


(99: 


which  can  be  solved  for  0  as 


9  =  sin  1 


Ap 


cos  a 


-e, 


(ioo: 


Once  again,  we  can  make  an  additional  verification  at  the 
boundary  between  sub-modes.  For  the  instantaneous  source 
voltage  vs  (0  =  Vm  sm(<yf3)  and  an  average  capacitor  voltage  as 
previously  solved  for  in  (70),  we  have  the  equality 


Vm  sin  Oh )  = 


APVm 


K 


cos  ab  . 


(ioi: 


By  substituting  cot3  =  0h+0L  into  (101),  we  have 


Vm  sin(04  +  9,  )  =  ^^cos ah , 
n 


(102: 


which  can  be  solved  for  0b  as 
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To  calculate  the  average  load  current,  we  again  refer  to 

(21).  By  Substituting  (108)  into  (21),  we  get 

r  pVma7 f  (n  ^  Ap(  9n  9,71  \  1J/  . 

1avg  =  ~ I  cos(0  +  0i)-cos(fi#)  + —  — +  cosa  d{cat)  (109) 

Ikcol  aJ+e  \_  n  'vl 80  180  )  J 

Because  of  the  complexity  of  (109),  it  is  again  broken  into 
three  parts.  Let 
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(110) 


a+0,. 


X1=cos(#  +  #i)  |  d(cot) , 


a+0. 


X,  =  X,  from  (82),  and 


a+0. 


Ap  cos  a  r'  On  0,n 
X,  =  — -  - +  -± - cot 

n  JAm  180  . 


3 


a  (cot) 


Starting  with  Xl  and  integrating  (110),  we  obtain 


X,  =  \jot  cos  (#  +  0L )] 


cot=a+0u 

cot=a+0L 


which,  when  evaluated,  becomes 


X, 


1  U80y 

which  can  be  rewritten  as 


{0U-0l)cos{0+0l)  , 


*i  = 


^2  n^ 
P  , 


COS 


(ff+0L) 


Because  X,  is  the  same  as  the  previous  X,  ,  its 
given  in  (91)  .  Finishing  with  X3  and  integrating 
obtain 


^3  = 


Ap  cos  a 


V  x  J 


On  0,n 


180  180 


cot- 


(| cot ) 


(ot=a+0,. 


cot=a+0, 


which,  when  evaluated,  becomes 

I  ^  On  0,n 5 

f  Ap  cos  a^ 


X3  = 


V  n  ) 


180  +  180 


(0u-0l) 


(  n  ^ 


V  1 80y 


1  \(  an  01Tn\  (  an  0,n 5 

-+  L 


180  180 


180  180 


By  simplifying  (116),  the  result  is 


(111) 

(112) 

(113) 

(114) 

result  is 
(111),  we 

(115) 


(116) 
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X3  =2  A 

By  inserting  the  results 
into  (109),  we  have 


^  On  an 

ym~m 


\ 

n 

P V 


cos  a  . 


from  (114),  (91), 


and  (117) 


/  =FLl 

avs  ncoL 


\PJ 


cos(6,  +  6'i)  +  y4sina  +  A 


On 


an  n 


180  180  p  J 


cos  a 


By  using  (97),  we  can  substitute  it  into  (118)  and 
for  Lcrit  as 


L.t  = 


R 


coAcosa 


f n ' 


\PJ 


cos  ( 6  +  0L )  +  A  sin  a  +  A 


On 


an  n 


180  180  p  J 


cos  a 


Although  0  is  a  function  of  a  as  given  by 
substituting  it  into  (119)  would  only  serve  to  make  i 
unwieldy,  so  we  leave  it  in  its  present  form. 

5 .  Summary  of  Calculations  for  Lcrit 


To  summarize,  we  have  the  following  equations  for 


L  = 

crit 


Rtana 


co 


f  x  nB ^ 
v  APj 


for  a>a. 


and 


L  = 


R 


coAcosa 


f  n^ 


\P) 


cos  {0  +  0L )  +  A  sin  a  +  A 


^  On  an  n^ 


180  180  p 


cos  a 


for  a  <ab ,  where 


ab  =  tan 


p  B 
n  A 


0  =  sin  1 


Ap 


cos  a 

v  n  j 


-o,, 


(117) 

back 

(118) 

solve 

(119) 

(100) , 
■  more 

T 

^ crit  # 

(120) 

(121) 

(122) 

(123) 
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(124) 


A  =  cos  0 L  , 

and 

B  =  sm0L.  (125) 

From  Figure  13,  critical  inductances  for  a  6-,  12-,  18-, 
24-,  and  96-pulse  rectifier  can  be  seen,  assuming  a  unity 
resistance  and  a  frequency  of  60  Hz.  As  before,  as  a 
increases  for  each  rectifier,  Lcrit  increases  at  an  increasing 
rate.  Another  important  conclusion  from  Figure  13  is  that 
as  pulse  count  increases,  Lcrit  decreases  at  a  decreasing 
rate.  This  is  most  obvious  when  comparing  a  6-pulse 
rectifier  to  a  12-pulse  rectifier,  where  there  is  a  large 
reduction  in  Lcrit ,  and  then  continuing  to  an  18-pulse 
rectifier,  where  the  percent  difference  in  Lcrit  is  less 
substantial.  To  clearly  illustrate  the  concept  of 
diminishing  returns,  the  jump  from  a  24-pulse  rectifier  to 
a  96-pulse  rectifier  was  purposely  plotted  to  show  just  how 
little  reduction  in  Lcrit  is  obtained  by  dramatically 
increasing  pulse  count.  Finally,  as  pulse  count  increases, 
ah  decreases,  making  the  equation  for  Lcrit  for  a  <ah  less 
relevant,  especially  for  high  pulse  count  rectifiers. 
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x  10' 


a  (degrees) 


Figure  15.  Critical  inductance  plot  for  controlled 
p-pulse  rectifiers  with  unity  resistance 
operating  at  60  Hz. 


Critical  inductance  values  for  the 

15  are  provided  in  Table  1,  starting  at 
0°  and  ending  at  80°,  in  10°  increments. 
90°  would  result  in  an  infinite  critical 
upper  bound  is  not  inclusive  of  90°. 


Lcrit  plots  of  Figure 
a  firing  angle  of 
,  A  firing  angle  of 
inductance,  so  the 
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Table  1 .  Critical  inductance  (Henry)  for  selected  firing 
angles  for  the  controlled  p-pulse  rectifiers  of 

Figure  15. 


firing  angle  alpha 

degrees) 

pulse-count 

0 

10 

20 

30 

40 

50 

60 

70 

80 

6 

2.51E-05 

4.35E-05 

8.99E-05 

1.43E-04 

2.07E-04 

2.94E-04 

4.28E-04 

6.79E-04 

1.40E-03 

12 

3.07E-06 

1.07E-05 

2.22E-05 

3.51E-05 

5.11E-05 

7.26E-05 

1.05E-04 

1.67E-04 

3.45E-04 

18 

9.07E-07 

4.76E-06 

9.82E-06 

1.56E-05 

2.26E-05 

3.22E-05 

4. 67  E- 05 

7.42E-05 

1.53E-04 

24 

3.82E-07 

2.67E-06 

5.52E-06 

8.76E-06 

1.27E-05 

1.81E-05 

2.63E-05 

4.17E-05 

8.60E-05 

96 

5.96E-09 

1.67E-07 

3.45E-07 

5.47E-07 

7.95E-07 

1.13E-06 

1.64E-06 

2.60E-06 

5.37E-06 

A  much  more  useful  table  would  be  one  that  shows  a 
relative  comparison  between  the  absolute  Lcrit  values  of 
different  controlled  p-pulse  rectifiers,  such  as  Table  2. 
In  Table  2,  the  6-pulse  controlled  rectifier  is  treated  as 
the  basis  to  which  the  other  p-pulse  controlled  rectifiers 
are  compared.  The  Lcrit  ratios  become  fixed  percentages  beyond 
the  boundary  angle  for  each  controlled  rectifier. 


Table  2.  Critical  inductance  as  a  percentage  of  6-pulse 
critical  inductance  for  selected  firing  angles  for 
the  controlled  p-pulse  rectifiers  of  Figure  15. 


firing  angle  alpha 

degrees) 

pulse-count 

0 

10 

20 

30 

40 

50 

60 

70 

80 

6 

100% 

100% 

100% 

100% 

100% 

100% 

100% 

100% 

100% 

12 

12% 

25% 

25% 

25% 

25% 

25% 

25% 

25% 

25% 

18 

4% 

11% 

11% 

11% 

11% 

11% 

11% 

11% 

11% 

24 

2% 

6% 

6% 

6% 

6% 

6% 

6% 

6% 

6% 

96 

<1% 

<1% 

<1% 

<1% 

<1% 

<1% 

<1% 

<1% 

<1% 
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IV.  DC  OUTPUT  VOLTAGE  RIPPLE 


The  discussion  of  critical  inductance  began  with  the 
assumption  that  the  instantaneous  capacitor  voltage  was 
approximately  equal  to  the  average  capacitor  voltage  (which 
was  equal  to  the  average  load  voltage) .  While  it  was  useful 
then  to  conceptually  think  of  the  capacitor  as  being 
infinitely  big  and  causing  vc(t )  to  be  constant,  we  know  this 
is  not  the  case.  The  load  voltage  will  have  some  amount  of 
ripple,  but  it  will  be  mitigated  by  the  capacitor  during 
its  charge  and  discharge  periods. 

A.  OUTPUT  VOLTAGE  RIPPLE  FOR  A  CONTROLLED  RECTIFIER  WITH 

A  RESISTIVE  LOAD 

While  the  concentration  of  this  thesis  is  not 
specifically  about  phase  controlled  rectifiers  with  purely 
resistive  loads,  it  is  a  good  starting  place  to  understand 
about  output  voltage  ripple.  The  circuit  of  Figure  16  is  an 
example  of  a  controlled  three-phase  (6-pulse)  rectifier 
with  only  a  resistive  load.  However,  let  us  consider  it  to 
be  a  controlled  p-pulse  rectifier  for  the  proceeding  output 
ripple  voltage  analysis. 
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Figure  16.  Circuit  diagram  of  a  controlled  three-phase 
rectifier  with  a  resistive  load.  From  [4] . 


From  [3]  and  [4],  for  a<0L,  the  instantaneous  load  voltage 
is  always  positive  or  zero,  and  its  average  load  voltage  is 
synonymous  to  (70)  and  is  given  by 

V  =  ApVm  cos  a  (12  6) 

n 

For  a>0L,  the  instantaneous  load  voltage  becomes  negative 
in  parts  of  the  cycle  and  is  expressed  as 

180° 

Vavg=^~  \  Vmsinrd/  r  (127) 

where  the  upper  integration  limit  is  replaced  with  180°  to 
limit  the  instantaneous  load  voltage  to  either  zero  or 
positive  values.  Integrating  and  evaluating  (127),  we  get 

,7flvg=V^[1  +  C0S(«  +  6'i)]  (12  8) 

In 
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From  [4],  the  root  mean  square  (RMS)  value  of  the  load 
voltage  is  obtained  by  applying  the  definition  of  RMS, 
given  by 


Vrms  \  I  j  j 


1  0 

Starting  with  the  first  case,  we  have 


:i29) 


V  = 

r  RMS 


pK 


2  a+Oy 


2  n 


J  sin 1  ydy 


( 130 ; 


(x+Qi 


By  applying  the  trigonometric  identity  skT0  =  [l-  cos(2#)]/2  to 
(130)  and  integrating,  we  get 


V  = 

'  RMS 


pK; 

8n 


1  r=a+ev 


{2r-sin(2y)]^ 


(131) 


which,  when  evaluated,  becomes 


pV  2 

m  J 

|4^_r 

8  TV 

L  P  ^ 

(132; 


By  applying  sm(a  +  b)  =  sin(o)cos(&)  +  cos(tf)sin(&) ,  we  have 


v  = 

v  RMS 


pVj 

87V 


—  +  sin  ( 2a )  cos  ( 26 L )  +  cos  ( 2a )  sin  ( 26 L  ) 

-  sin  ( 2a )  cos  ( 26^ )  -  cos  ( 2a )  sin  ( 26^ ) 


(133) 


Next,  let  E  =  cos(26l),  F  =  sin(2dL)  ,  G  =  cos(2 0V)  ,  and  F[  =  sm{26u)  ; 
therefore,  H  —  —F  and  G  =  E.  Then,  (133)  reduces  to 


_  Vm  Fp  cos  (2a) 
Vrms  2n 


(134; 


Continuing  with  the  second  case,  we  have 
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(135) 


v  = 

'  RMS 


pK 


2  180° 


2  n 


J  sin^  ydy 


a+6. 


We  again  apply  sin2  6  —  [l- cos(20)]/ 2  to  this  result 
integrate  to  obtain 


V  = 

r  RMS 


PVm 

8n 


-[2y-sin(2y)] 


y=180° 

y=a+6L 


which,  when  evaluated,  becomes 


v 

r  RMS 


\pVm 


8  n 


Ct7Z  2n  .  /  „  \ 

k - +  —  +  sin  ( 2  a  +  20, 

90  p  y  ’ 


By  again  applying  sin(a  +  b)  =  sin(tf)cos(6)  +  cos(a)sin(&) ,  we  have 


V  - 

v  RMS 


PLL 

8n 


„  an  2  n 

2  n - n-\ - 

90  n 


+  sin  (2«)  cos  [26 L )  +  cos  (2a)  sin  [20 L ) 


Next,  let  E  =  cos(26l)  and  F  =  sin(2#t).  By  substituting 
into  (138)  and  rearranging,  we  get 


_Vm  2  +  p  ap  Epsm{2a)  Fpcos(2a) 
4  360 +  4 n  +  4n 


Now  we  can  compute  the  RMS  ripple  voltage  Epms-rippie  for 
case  as 


v  =  [v1  -V2 

v  RMS— ripple  \]  v  RMS  v  avg 


For  the  first  case  ( a<0L ),  we  substitute  (134)  and 
into  (140)  to  obtain 


and 

(136) 

(137) 

(138) 

these 

(139) 

each 

(140) 

(126) 


46 


Computing  the  peak-to-peak  ripple  voltage  Vpp_ri  le  and  the 
peak-to-peak  ripple  voltage  as  a  percentage  of  the  average 
output  voltage  %Vpp/avg  are  slightly  more  involved  because 

there  are  actually  three  separate  regions  we  must  evaluate. 
The  first  two  regions  are  actually  sub-regions  of  the 
aforementioned  first  case,  while  the  final  region 
correlates  with  the  second  case  from  earlier.  For  the  first 

region,  given  by  a  <  {0U  —  0L^/ 2 ,  the  maximum  voltage  is  Vm , 

and  the  minimum  voltage  is  sin  (a +  6^)  .  Therefore,  Vpp-ripplei  is 

their  difference, 

Vpp-nppie,  =  K,  [l  -  sin  ( a  +  0U )]  ,  (14  3) 

which  can  be  written  in  terms  of  p  as 


(  180°^) 

Vpp-ripple=K  1-sin  a+  90°  + -  (144) 

v  P  )_  ■ 

The  peak-to-peak  ripple  voltage  as  a  percentage  of  the 
average  output  voltage  is  given  by 

( V  N 

%Vpp/avg  =  100%.  (145) 
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By  substituting  (144)  and  (126)  into  (145),  we  get 


1-sin  a +  90°  +  - 


100% 


(146; 


The  second  region  is  given  by  [0u—6^l2<a<6L,  in  which  the 
maximum  voltage  is  Vmsi n(a  +  0L^,  while  the  minimum  voltage  is 
Vmsm(a  +  0u')  .  Therefore,  Vpp_ripple j  for  this  second  interval  is 


pp-ripple 2 


■■  Fm[sin(a  +  0i)-sin(a  +  0c/)]  , 


which  can  be  written  in  terms  of  p  as 


i  1 80°  i  (  1 80° 

Vpp.ripplei=Vm  Sin  a  +  90° - sin  a +  90°  + - 

V  p  )  V  P 


By  substituting  (148)  and  (126)  into  (145),  we  get  the 
peak-to-peak  ripple  voltage  as  a  percentage  of  the  average 
output  voltage 


I  1 80°  i  (  1 80° 

sin  a  +  90° - sin  a  +  90°  + - 


100% 


(149; 


Finally,  the  third  region  is  given  by  a>dL.  In  this  region, 
the  maximum  voltage  is  Vmsin(a  +  ,  and  the  minimum  voltage 
is  zero,  making  Vpp_ripple  equal  to  the  maximum  voltage  or 


V pp-ripple,  =VmSin(a  +  0L)r 


(iso; 


which  can  be  written  in  terms  of  p  as 
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(151) 


f 


V  =  V 

pp  -ripple^  m 


sin 


a  +  90° 


180°^ 

P  y  • 


By  substituting  (151)  and  (128)  into  (145),  we  get  the 
peak-to-peak  ripple  voltage  as  a  percentage  of  the  average 
output  voltage 


%V, 


pp/avg3 


f 


sin 


a  +  90°- 


180c 

P  J 


^-[l  +  cos(a  +  &L)\ 
2.71 


100% 


(152) 


The  results  of  both  cases  for  a  6-pulse  rectifier  with 
Vm  egual  to  1.0  V  and  a  purely  resistive  load  are 
summarized  in  Figures  17  and  18.  Initially,  as  a  increases, 
the  difference  between  the  RMS  load  voltage  and  the  average 
load  voltage  increases,  making  the  sguare  root  of  the 
difference  of  their  squares  grow,  causing  the  RMS  load 
voltage  ripple  to  steadily  increase.  However,  as  a 
continues  beyond  60°,  the  square  root  of  the  difference  of 
their  squares  begins  to  decrease,  causing  the  RMS  load 
voltage  ripple  to  peak  and  also  begin  to  decrease.  The 
peak-to-peak  ripple  load  voltage  rapidly  increases  until 
reaching  a  subtle  inflection  point  at  a  =  30°  and  then 
increases  less  rapidly  until  reaching  a  maximum  at  a-  60°, 
where  it  begins  to  rapidly  decrease. 
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0  20  40  60  80  100  120 

a  (degrees) 

Figure  17.  Average  voltage,  RMS  voltage,  RMS  ripple 

voltage,  and  peak-to-peak  ripple  voltage  for  a  6-pulse 
rectifier  circuit  with  a  unity  Vm  and  a  purely 

resistive  load. 

The  peak-to-peak  ripple  voltage  as  a  percentage  of  the 
average  output  voltage  increases  as  a  increases  to  60°  (the 
transition  at  30°  is  virtually  undetectable)  because  Vavg  is 
fairly  high  and  decreasing  while  V  _rt  le  is  low  but  growing 
quickly.  Beyond  60°,  %Vpp/avg  increases  at  a  tremendous  rate 
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This 


due  to  V  decreasing  at  a  faster  rate  than  Vpp_ripple 

illustrates  another  reason  why  circuit  designers  prefer  to 
keep  a  fairly  small. 


Figure  18.  Peak-to-peak  ripple  voltage  as  a 
percentage  of  the  average  output  voltage  for  a  6-pulse 
rectifier  circuit  with  a  unity  Vm  and  a  purely 

resistive  load. 

In  order  to  demonstrate  how  increasing  pulse  count 
improves  the  peak-to-peak  ripple  voltage  as  a  percentage  of 
the  average  output  voltage  for  these  types  of  rectifiers, 
%Vpp,aVg  has  been  plotted  for  several  p-pulse  rectifier 


circuits  in  Figure  19.  The  plot  exhibits  very  similar 
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behavior  to  Figure  15  in  that  increasing  pulse  count  yields 
better  results  (with  rapidly  diminishing  returns) ,  while 
increasing  a  worsens  the  results.  While  not  shown  in  Figure 
19,  further  increasing  a  beyond  0L  (or  much  beyond  30°,  for 
that  matter)  is  purely  academic  as  %Vpp/  very  rapidly 

reaches  exorbitant  amounts.  This  also  illustrates  why  a 
rectifier  with  a  purely  resistive  load  is  an  unlikely 
candidate  for  a  rectifier  which  necessitates  a  low  %V  .  as 

it  will  require  a  high  pulse  count  (increased  hardware 
complexity,  cost,  and  space) . 


Figure  19.  Peak-to-peak  ripple  voltage  as  a  percentage 
of  the  average  output  voltage  for  p-pulse 
rectifier  circuits  with  unity  Vm  and 
purely  resistive  loads. 
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To  help  in  understanding  some  of  the  items  associated 
with  angles  for  the  p-pulse  rectifiers  with  purely 
resistive  loads,  a  quick  reference  summary  is  provided  in 
Table  3.  For  instance,  for  a  12-pulse  rectifier  with  a 
purely  resistive  load,  its  first  pulse  starts  at  75°  and 
ends  at  105°,  making  its  pulse  width  30°  (their 
difference) ;  its  three  regions  of  varying  peak-to-peak 
voltage  (and,  hence,  differing  %Vpp/avg )  are  the  intervals  of 
a  from  0°  to  15°,  15°  to  75°,  and  75°  to  105°. 

Table  3.  Angular  characteristics  of  some  p-pulse 

controlled  rectifiers  with  purely  resistive  loads. 


0L 

6u 

pulse-width 

Ot  limit  1 

Otlimit  2 

Supper  limit 

pulse-count 

(degrees) 

(degrees) 

(degrees) 

(degrees) 

(degrees) 

(degrees) 

6 

60 

120 

60 

30 

60 

120 

12 

75 

105 

30 

15 

75 

105 

18 

80 

100 

20 

10 

80 

100 

24 

82.5 

97.5 

15 

7.5 

82.5 

97.5 

96 

88.125 

91.875 

3.75 

1.875 

88.125 

91.875 

Just  as  Table  1  was  developed  to  provide  Lcrit  for 
selected  values  of  a  ,  Table  4  does  the  same  for  %V ,  for 

the  waveforms  of  Figure  19.  Firing  angles  of  120°  for  the 
6-pulse  controlled  rectifier  and  100°  for  the  18-pulse 
controlled  rectifier  would  result  in  infinite  %Vpp/avg  because 

those  values  of  a  are  the  upper  limits  (and,  therefore,  not 
inclusive)  for  the  two  rectifiers,  in  accordance  with  Table 

3. 
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Table  4 . 


Peak-to-peak  ripple  voltage  as  a  percentage  of 
the  average  output  voltage  for  selected  firing 
angles  for  the  controlled  p-pulse  rectifier 
waveforms  of  Figure  19. 


firing  angle  alp 

la  (degrees) 

pulse-count 

0 

10 

20 

30 

40 

50 

60 

70 

80 

90 

100 

110 

6 

14% 

25% 

40% 

60% 

88% 

125% 

181% 

225% 

288% 

391% 

594% 

1197% 

12 

3% 

10% 

19% 

30% 

44% 

62% 

91% 

144% 

236% 

398% 

1199% 

18 

2% 

6% 

13% 

20% 

29% 

42% 

60% 

96% 

198% 

399% 

OO 

24 

<1% 

5% 

10% 

15% 

22% 

31% 

45% 

72% 

148% 

399% 

96 

<1% 

1% 

2% 

4% 

5% 

8% 

11% 

18% 

37% 

400% 

While  the  data  presented  by  Table  4  is  somewhat  useful 
for  quickly  finding  the  value  of  %V  .  for  a  given  value  of 

a  for  a  particular  p-pulse  controlled  rectifier  with  a 
resistive  load,  another  way  to  look  at  the  data  is  by 
assigning  the  %Vpplavg  of  the  6-pulse  controlled  rectifier  as 

the  basis  to  which  the  others  are  compared.  The  result  is 
Table  5.  For  realistic  values  of  a,  higher  pulse  count 
controlled  rectifiers  clearly  possess  much  less  %V  ,  than 
lower  pulse  count  variants. 
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Table  5 . 


Peak-to-peak  ripple  voltage  as  a  percentage  of 
the  average  output  voltage  as  a  percentage  of  6- 
pulse  %V  .  for  selected  firing  angles  for  the 

controlled  p-pulse  rectifier  waveforms  of  Figure 

19. 


iring  angle  a 

pha  (c 

egrees) 

pulse-count 

0 

10 

20 

30 

40 

50 

60 

70 

80 

90 

100 

110 

6 

100% 

100% 

100% 

100% 

100% 

100% 

100% 

100% 

100% 

100% 

100% 

100% 

12 

25% 

39% 

48% 

50% 

50% 

50% 

50% 

64% 

82% 

102% 

202% 

18 

11% 

25% 

32% 

33% 

33% 

33% 

33% 

43% 

69% 

102% 

oo 

24 

6% 

19% 

24% 

25% 

25% 

25% 

25% 

32% 

52% 

102% 

96 

<1% 

5% 

6% 

6% 

6% 

6% 

6% 

8% 

13% 

102% 

B.  OUTPUT  VOLTAGE  RIPPLE  FOR  A  CONTROLLED  RECTIFIER  WITH 
A  TWO- POLE  LC  OUTPUT  FILTER 

Turning  back  to  the  circuit  of  Figure  9,  we  now  desire 

to  calculate  the  output  ripple  for  the  type  of  circuit  we 

are  primarily  interested  in  analyzing.  It  turns  out  that 
the  voltage  ripple  calculations  for  a  controlled  rectifier 
with  a  two-pole  LC  output  filter  are  actually  somewhat 

simpler  than  those  for  the  controlled  rectifier  with  a 
purely  resistive  load  because  we  are  going  to  assume  the 
inductor  is  always  large  enough  to  maintain  a  continuous 
current;  therefore,  there  is  only  one  case  to  evaluate,  not 
two.  Additionally,  some  of  the  calculations  are  the  same  as 
for  the  purely  resistive  case,  so  we  need  not  repeat  our 
previous  work.  To  begin  from  [4],  the  average  capacitor 
voltage  is  given  by  (70) ,  and  the  RMS  load  voltage  is  given 
by  (134)  .  Then,  the  RMS  ripple  voltage  is  given  by 

V  =  Iv2  -V2  (153) 

r  RMS-ripple  Sj  r  RMS  r  C(avg )  V  -L  ^  / 
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Here  we  will  make  a  small  notation  change,  letting  VDC  be 
equal  to  V  and  V AC  be  equal  to  k  to  better  illustrate 

the  relationship  between  RMS,  DC,  and  AC  quantities. 
Additionally,  it  is  important  to  understand  that  VAC  is 
applied  across  the  entire  load  (which  includes  the 
resistor,  capacitor,  and  the  inductor).  Therefore,  (153) 
becomes 

V  -  V2  (154) 

r  AC  V  RMS  r  DC  '  ±  U  ^  > 

and  (141)  becomes 


vAC  =  vj- 


1  + 


Fp  cos  (2a) 


2  n 


A2p2  2 
- cos  a 


(155) 


n 


The  impedance  of  the  capacitor  is 

Zc  s  — - — 

jpaC 


(156) 


and  the  impedance  of  the  inductor  is 

ZL  =  jpcoL  ,  (157) 

where  these  are  considered  approximations  because  Zc  and  ZL 
are  composed  of  the  frequencies  pco ,  2pco ,  3pco ,  etc.; 
however,  the  dominant  frequency  is  pco  .  Next,  the  impedance 
of  the  capacitor  in  parallel  with  the  resistor  is 


Z  - 

^R\\C  - 


R 


1  +  jpcoRC 


(158) 


which  can  be  written  in  terms  of  real  and  imaginary  parts 
as 


«||C 


R 

1  +  {pcoRC)2 


pcoRC 

’  l  +  (pcoRC)2  ■ 


(159) 
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The  impedance  of  ZL  in  series  with  ZR„C  is 


R 


Zt0,al  1  +  ( pcoRC) 
which  can  be  rewritten  as 

R 


-  +  j 


pcoL  - 


pcoR  C 
1  +  (pooRC)2 


7  - 

total 


■  1 


pcoL  +  p3a>3R2LC2  -  pcoR2C 


1  +  [pcoRC) 

The  magnitude  of  Zlolal  becomes 


1  +  [pcoRC)2 


7  - 

total  —  ■ 


r  r 

2 

_l_ 

pcoL  +  p2 a? R2LC2  -pcoR2C 

[1  +  (pcoRC)2 

i 

\  +  {pcoRC )2 

We  can  now  solve  for  Vac\  from  Ohm' s  Law  as 


,  ,  v 

\j  ~  v  AC 

AC  \r-y 


J  total 


By  substituting  (155)  and  (162)  into  (163),  we  get 


I  AC  = 


r  r 

2 

_|_ 

pcoL  +  p3o3R2LC2  -  pooR2C 

[1  +  (pcoRC)2 

i 

1  +  (pcoRC)2 

Therefore,  the  capacitor  RMS  voltage  VC(RMS)  is  given  by 


V  -I  Z 

r  C(RMS)  ~  1  AC  ^C 


and  the  capacitor  peak-to-peak  ripple  voltage  is 


V. 


C  ( pp-fipple) 


(is! 2) 


V, 


C(RMS) 


(160) 


(161) 


(162) 


(163) 


(164) 


(165) 


(166) 
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Here,  it  is  important  to  note  another  approximation  that 


has  taken  place.  We  have  made  the  assumption  that  Vt 


C(RMS) 


is  a 


perfect  sinusoid  such  that  its  peak  value  is  simply  equal 

V2F, 


to  '  “  C(RMS) 

voltage  as 

%VC(pP,avg)  IS 


Finally,  the  capacitor  peak-to-peak  ripple 
a  percentage  of  the  average  output  voltage 


%V, 


C(pplavg) 


V 

C(pp-ripple) 

V  ^C(avg)  ) 


100%  , 


(i67: 


where  VC(avg)  is  given  by  (70)  and  VC(pp_ripple)  is  given  by  (166)  . 

One  of  the  things  that  would  be  interesting  to 
investigate  is  how  %V„,  ,  behaves  as  the  individual 

3  C(pplavg) 

circuit  elements  (inductor,  capacitor,  and  resistor)  are 
changed  one  at  a  time.  While  the  behavior  of  0//°^c(pp/avg)  for 
every  circuit  is  unique  depending  on  its  peak  line-to-line 
input  voltage  Vm  ,  pulse  count,  and  the  values  of  the 
inductor,  capacitor,  and  resistor,  we  will  examine  one  case 
with  the  aid  of  Figure  20.  The  particular  circuit  that 
generated  the  waveforms  of  Figure  20  is  a  6-pulse 
controlled  rectifier  (in  the  configuration  given  by  Figure 
9)  with  a  Vm  of  100  V,  a  25  mH  inductor,  a  50  |iF  capacitor, 

and  a  10Q  resistor.  Through  some  initial  trial  and  error 
with  different  resistance  values,  it  was  discovered  that 
%VC(pp/avg)  was  very  insensitive  to  changes  in  the  resistance, 

so  it  was  excluded  from  further  sensitivity  analysis.  The 
initial  values  of  the  inductor  and  capacitor  are  as 
previously  mentioned  and  are  designated  as  L0  and  C0  in  the 
legend  of  Figure  20. 
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Figure  20.  Capacitor  peak-to-peak  ripple  voltage 

as  a  percentage  of  the  average  output  voltage  for  a  6- 
pulse  controlled  rectifier  with  a  two-pole  LC  output 

filter . 

Somewhat  arbitrarily,  it  was  decided  to  increase  and 
decrease  L0  and  C0  by  fifty  percent  to  observe  the  effect  on 
%F„,  ,  v  (care  had  to  be  taken  to  ensure  the  circuit  was 

C(pplavg) 

still  in  CCM  for  the  reduced  inductor  value)  .  The  larger 
inductor  and  capacitor  are  designated  as  LH  and  CH  in  the 
legend  of  Figure  20,  while  the  reduced  components  are  Ll 
and  CL  ,  respectively.  Interestingly,  the  combination  of  a 
reduced  inductor  and  initial  capacitor  (worst  0//°^c(pp/mg) ) 

produced  close  results  to  the  initial  inductor  and  reduced 
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capacitor  pair,  while  the  larger  inductor  and  initial 
capacitor  produced  nearly  identical  results  to  the 
combination  of  the  initial  inductor  and  bigger  capacitor— 
both  of  which  could  be  said  to  have  the  best  %Vrl  ,  ,  . 

In  order  to  see  how  increasing  pulse  count  affects 
0/°Vc(pp/avg) '  °^c{ppiavg)  was  plotted  in  Figure  21  for  a  6-,  12- , 

18-,  24-,  and  96-pulse  controlled  rectifier  with  a  two-pole 

LC  output  filter  using  the  same  initial  component  values 
and  Vm  as  before. 

m 


a  (degrees) 


Figure  21.  Capacitor  peak-to-peak  ripple  voltage 

as  a  percentage  of  the  average  output  voltage  for  a 
6-,  12-,  18-,  24-,  and  96-pulse  controlled  rectifier 
with  a  two-pole  LC  output  filter. 
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Much  like  the  plots  of  Figures  15  and  19,  tremendous 
improvement  is  noted  in  increasing  from  a  6-pulse  to  a  12- 
pulse  controlled  circuit,  but  further  increasing  pulse 
count  results  in  less  of  an  improvement. 

The  values  of  %Vr(  ,  ■>  for  selected  values  of  a  for  the 

C(pp/avg) 

waveforms  of  Figure  21  are  collected  in  tabular  form  in 
Table  6.  The  capacitor  peak-to-peak  ripple  voltage  as  a 
percentage  of  the  average  output  voltage  approaches 
infinity  at  the  firing  angle  limit  of  90°  because  FC(  }  is 

zero  at  90°;  hence,  the  firing  angle  limit  is  exclusive  of 
90°  . 


Table  6.  Capacitor  peak-to-peak  ripple  voltage  as  a 
percentage  of  the  average  output  voltage  for 
selected  firing  angles  for  the  controlled  p-pulse 
rectifier  waveforms  of  Figure  21. 


firing  angle  alpha  (degrees) 

pulse-count 

0 

10 

20 

30 

40 

50 

60 

70 

80 

6 

2% 

3% 

6% 

9% 

13% 

18% 

26% 

41% 

84% 

12 

<1% 

<1% 

<1% 

1% 

1% 

2% 

3% 

5% 

10% 

18 

<1% 

<1% 

<1% 

<1% 

<1% 

<1% 

<1% 

1% 

3% 

24 

<1% 

<1% 

<1% 

<1% 

<1% 

<1% 

<1% 

<1% 

1% 

96 

<1% 

<1% 

<1% 

<1% 

<1% 

<1% 

<1% 

<1% 

<1% 

As  stated  previously,  every  circuit  of  this  type  is 
unique  and  will  produce  differing  %Vc^pp/avgj  depending  on  the 

parameters  of  the  controlled  rectifier  with  a  two-pole  LC 
output  filter;  therefore.  Table  6  has  limited  utility. 
However,  if  we  employ  the  technique  we  have  previously  used 
and  treat  %VC(pp/avg)  f or  the  6-pulse  controlled  rectifier  as 
the  basis  to  which  the  other  rectifiers  are  compared,  we 
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generate  Table  7.  As  opposed  to  Table  6,  which  deals  in 
absolutes,  the  relationships  in  Table  7  are  always  valid  no 
matter  what  the  controlled  rectifier  circuit  parameters  and 
value  of  Vm  are  (so  long  as  the  circuit  remains  in  CCM) . 


Table  7 .  Capacitor  peak-to-peak  ripple  voltage  as  a 
percentage  of  the  average  output  voltage  as  a 
percentage  of  6-pulse  %VC(pp!avg)  for  selected  firing 

angles  for  the  controlled  p-pulse  rectifier 
waveforms  of  Figure  21. 
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100% 

100% 

100% 

100% 

100% 

100% 

100% 

100% 

100% 

12 

5.8% 

9.9% 

11.1% 

11.4% 

11.6% 

11.6% 

11.6% 

11.7% 

11.7% 

18 

1.1% 

2.8% 

3.2% 

3.3% 

3.4% 

3.4% 

3.4% 

3.4% 

3.4% 

24 

<1% 

1.1% 

1.3% 

1.4% 

1.4% 

1.4% 

1.4% 

1.4% 

1.4% 

96 

<1% 

<1% 

<1% 

<1% 

<1% 

<1% 

<1% 

<1% 

<1% 
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V.  HARMONICS 


A.  BACKGROUND 

Harmonics  have  existed  on  power  systems  from  the  time 
of  the  very  first  generators.  At  that  time,  the  harmonic 
components  were  so  small  that  their  effects  on  systems  were 
negligible.  This  was  due  to  the  lack  of  non-linear  loads 
before  the  1960s  [12]  .  However,  in  recent  years,  there  has 
been  a  significant  increase  in  the  installation  and  use  of 
power  electronics  equipment  onboard  ships  and  at  offshore 
installations.  The  operation  of  this  equipment  can,  and  has 
in  many  cases,  significantly  degraded  the  ship  or  offshore 
installation  electrical  power  quality  to  such  an  extent 
that  corrective  measures  have  to  be  implemented  in  order  to 
reduce  the  resultant  adverse  effects  on  the  electrical 
plant  and  equipment.  The  quality  and  security  of  voltage 
supplies  are  important  to  the  safety  of  any  vessel  and  its 
crew  and  to  the  protection  of  the  marine  environment.  Any 
failure  or  malfunction  of  equipment  such  as  propulsion  or 
navigation  systems  can  result  in  an  accident  at  sea  or 
close  ashore  with  dire  consequences.  What  we  are  talking 
about  is  harmonic  distortion  of  voltage  supplies  caused  by 
the  operation  of  electronic  devices  which  draw  non-linear 
(i.e.,  non-sinusoidal  in  nature)  currents  from  the  voltage 
supplies.  These  very  same  pieces  of  non-linear  equipment 
can  also  be  adversely  affected  by  harmonic  currents  and  the 
subsequent  voltage  distortion  they  produce,  as  can  the 
majority  of  linear  equipment  (particularly  generators,  AC 
motors  and  transformers)  [13] . 
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B.  COMMON  PROBLEMS  ASSOCIATED  WITH  HARMONICS 

Power  converters  are  considered  the  primary  source  of 
undesired  harmonics.  These  AC/DC  power  converters  are 
extensively  used  in  various  applications  like  power 
supplies,  DC  motor  drives,  front-end  converters  in 
adjustable-speed  AC  drives,  high  voltage  DC  transmission, 
switched-mode  power  systems,  fluorescent  lights  (with 
electronic  ballasts) ,  utility  interfaces  with  non- 
conventional  energy  sources,  in-process  technology  like 
welding,  power  supplies  for  telecommunications  systems, 
aerospace,  and  military  environments.  Traditionally,  AC/DC 
power  conversion  has  been  dominated  by  diode  or  phase 
controlled  rectifiers  which  act  as  non-linear  loads  on 
power  systems  and  draw  input  currents  that  are  rich  in 
harmonics  and  have  poor  supply  power  factor,  thereby 
creating  the  power  quality  problem  for  the  distribution 
network  and  for  other  electrical  systems  in  the  same 
vicinity  of  the  rectifier.  Some  common  problems  associated 
with  the  production  of  these  harmonics  include: 

•  Large  reactive  power  drawn  by  rectifiers  from  the 
power  system  which  requires  that  the  distribution 
equipment  handle  large  power,  thus,  increasing 
its  volt-ampere  ratings; 

•  Voltage  drops  at  the  buses; 

•  Higher  input  current  harmonics  resulting  in  the 
distorted  line  current  which  tends  to  distort  the 
line  voltage  waveform.  This  often  creates 
problems  in  the  reliable  operation  of  sensitive 
equipment  operating  on  the  same  bus; 

•  Increased  losses  in  the  equipment  such  as 
transformers  and  motors  connected  to  the  utility; 

•  Electromagnetic  interference  with  any  nearby 
communications  circuits; 
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in 


•  Excessive  neutral  current,  resulting 
overheated  neutrals.  The  currents  of  triplen 
harmonics,  especially  odd  harmonics  (3rd,  9th 
15th,  etc.)  are  actually  additive  in  the  neutral 
of  three-phase,  wye-connected  circuits; 

•  Measuring  equipment  and  meters  reading 
incorrectly,  including  induction  disc-type  Watt- 
hour  meters  and  averaging  type  current  meters; 

•  Blown  fuses  on  power  factor  correction  capacitors 
due  to  high  voltages  and  currents  from  resonance 
with  line  impedance  and  capacitor  bank  failures; 

•  Mal-operation  of  equipment  such  as  computers, 
telephone  systems,  and  electronic  controllers; 

•  Nuisance  operation  of  protective  devices 
including  false  tripping  of  relays  and  failure  of 
uninterruptible  power  supplies  to  transfer 
properly,  especially  if  the  controls  incorporate 
zero-crossing  sensing  circuits;  and 

•  Damaging  dielectric  heating  in  cables  [14] . 

C.  VISUAL  REPRESENTATIONS  OF  HARMONICS 

Before  delving  directly  into  the  mathematics  and 
equations  dealing  with  harmonics,  it  is  useful  to  go  a  bit 
deeper  into  what  is  meant  by  a  load  being  either  linear  or 
non-linear.  A  linear  load  will  draw  current  that  is 
proportional  to  the  applied  voltage,  as  in  Figure  22  [13] . 


Voltage  and  current  waveforms  for  a 
linear  load.  From  [13] . 
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Figure  22 . 


It  stands  to  reason  then,  that  a  non-linear  load  will  draw 
a  current  that  does  not  flow  in  proportion  to  the  applied 
voltage  [15]  .  One  such  example  of  a  non-linear  load  is  the 
load  current  profile  for  a  three-phase  AC  drive  shown  in 
Figure  23. 


Figure  23.  Input  line  to  neutral  voltage  (in  black) 

and  input  current  (in  green)  of  phase  A  for  a  three- 
phase  AC  drive.  From  [15] . 

Around  the  1830s,  a  French  mathematician  named  Fourier 
discovered  that  a  distorted  waveform— like  the  current 
waveform  shown  in  Figure  23— can  be  represented  as  a  series 
of  sine  waves,  each  an  integer  number  multiple  of  the 
fundamental  frequency  and  each  with  a  specific  magnitude. 
For  example,  for  a  60  Hz  system,  the  fundamental  frequency 
is  60  Hz  and  the  5th  harmonic  would  be  five  times  60  Hz  or 
300  Hz  [15]  .  When  all  harmonics  are  added  to  the 
fundamental,  a  waveform  known  as  a  complex  wave  is  formed. 
A  simple  example  of  a  complex  wave  consisting  of  the 
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fundamental  and  3rd  harmonic  is  shown  in  Figure  24.  The 
complex  wave  is  clearly  the  algebraic  summation  of  the 
other  two  waves  [13] . 

60  Hz  Fundamental 


Figure  24.  Construction  of  a  complex  wave  from 

the  fundamental  and  3rd  harmonic.  From  [13] . 

The  collective  sum  of  the  fundamental  and  each 
harmonic  is  called  a  Fourier  series.  This  series  is  often 
viewed  as  a  spectrum  analysis  where  the  fundamental 
frequency  and  each  harmonic  component  are  displayed 
graphically  in  a  bar  chart  format  as  shown  in  Figure  25. 
The  total  current  is  the  square  root  of  the  sum  of  the 
squares  of  the  individual  harmonic  components.  One  general 
thing  we  can  say  regarding  the  harmonic  content  of 
something  is  by  simply  looking  at  the  wave  shape.  The 
more  it  looks  like  a  sine  wave,  the  lower  the  harmonic 
content  [ 15 ] . 
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Figure  25.  Harmonic  spectrum  analysis.  From  [15] . 


D.  FOURIER  SERIES  AND  TOTAL  HARMONIC  DISTORTION 

From  [16],  in  general,  any  non-sinusoidal  waveform  f(t) 
repeating  with  an  angular  frequency  co  can  be  expressed  as 

\ah  cos(hcot)  +  bhsm(hcot)}  ,  (168) 

h= i  2  h=i 


where  h  is  the  harmonic  number  and  F0  =  a0/2  is  the  average 
value.  Additionally,  we  have 


and 


1  2r 

ah  =—  /  (t)  cos  (hcot)d  [cot) 
n  i 


(169) 


K 


1  2r 

—  f(t)sin[hojt)d [cot)  , 
n  o 


(170) 


in  which  h  =  0,...,oo  .  From  (168)  and  (169),  we  can  solve  for 
the  average  value  F0  (taking  note  that  co  =  2kIT)  as 


]_ 

2 


Fn  =  —an  =  -3-2J  /(*)</ (of)  =  ^ }/(0* 


2;r 


(171) 


Oftentimes,  AC  waveforms  have  a  zero  average  value. 
Moreover,  by  the  use  of  waveform  symmetry,  it  is  often 
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possible  to  simplify  the  calculations  of  ah  and  bh  in  (169) 
and  (170)  .  A  list  of  the  types  of  symmetry,  requirements, 
and  formulas  for  ah  and  bh  is  provided  in  Table  8. 


Table  8.  Use  of  symmetry  in  Fourier  analysis.  From  [16] . 


Symmetry 

Condition  Required 

ah  and  bh 

Even 

1 

II 

2  fit 

bh  =  0 

ah  =  —  fit)cos(h(xit)  diwt) 

it  J0 

Odd 

f{~t)  =  -fit) 

2  fir 

o 

II 

<3 

bh  =  —  fU)sin(hwt)  di(ot) 

tt  Jo 

Half-wave 

fit)  =  -fit  +  \T) 

ah  =  bh 

=  0  for  even  h 

2 

f 11 

ah  =  ~ 

fit)  cos(/kj>7)  dibit)  for  odd  h 

17 

>0 

2 

f 11 

K  =  - 

fit)sinih(bt)  dibit)  for  odd  h 

17 

>0 

Even 

Even  and  half-wave 

<3- 

> 

II 

O 

for  all  h 

quarter-wave 

f 

—  f{t)co&ih(nt)  dibit)  for  odd  h 

tt  J0 

0  for  even  h 

Odd 

Odd  and  half-wave 

& 

a- 

II 

O 

for  all  h 

quarter-wave 

b»  = 

—  j  fit)s'm(huit)  dibit)  for  odd  h 

1 

0  for  even  h 

For  simplicity,  let  us  now  assume  that  we  have  a 
purely  sinusoidal  utility  input  voltage  at  the  fundamental 

frequency  (with  cox=co  and  fx  =  f  )  as  vs  (f)  =  y/2Vs  sin^f  .  Then, 

from  [16],  the  input  current  in  steady  state  is  the  sum  of 
its  Fourier  (harmonic)  components,  given  by 

*,(0=*,i(0+2X*(0/  (i72) 

h^\ 
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where  we  have  assumed  there  is  no  DC  component  in  is  ,  isl  is 
the  fundamental  (line  frequency  fl)  component,  and  ish  is  the 
component  at  the  h  harmonic  frequency  fh{=hfi)  •  The  amount  of 

distortion  in  the  current  waveform  is  quantified  by  means 
of  an  index  termed  the  THD.  The  distortion  part  of  the 
current  idis  is  given  by 

hiS (0  =  4 (0  -  4i (0  =  Z 4 h(t)  ■  ( 1 7  3 ) 

h±\ 

In  terms  of  RMS  values,  we  have 

-  <i74> 

where  the  RMS  value  of  is(t)  has  been  determined  by 


/  =  - 
\  T 


M[4(0  fdt 


(175) 


Finally,  the  THD  of  the  current  is  defined  as 


%THD ,  =  100% 


f  1  ^ 
1  dis 

V  7  S'!  J 


=  100% 


■1^: 


:  i7  6) 


E.  REDUCTION  OF  HARMONICS  BY  INCREASING  PULSE  NUMBER 

One  of  the  most  important  benefits  gained  by 
increasing  the  pulse  count  of  a  controlled  rectifier  is  a 
reduction  in  harmonic  content.  Without  diving  into  the 
detailed  mathematics  about  why  this  is  the  case,  suffice  it 
to  say  that  when  multiple  properly  phase  shifted  converters 
are  connected  in  series  or  parallel,  harmonic  cancellation 
occurs  such  that  the  line  currents  drawn  by  a  p-pulse 
converter  system  have  a  harmonic  content  given  by 
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h  =  kp±  1 , 


(177) 


where  k  =  1,2,3...  [6].  Additionally,  from  [16],  the  RMS  values 
of  the  harmonic  components  are  inversely  proportional  to 
their  harmonic  order  as 

(178) 

From  (177)  and  (178),  it  is  easy  to  understand  that  higher 
pulse  count  systems  will  have  smaller  harmonics  that  are 
further  away  from  the  fundamental  component. 

It  is  natural  to  consider  a  6-pulse  controlled 
rectifier  system  as  the  basic  building  block  from  which 
higher  pulse  count  systems  are  constructed.  In  order  to 
achieve  the  desired  harmonic  elimination  using  stacked  6- 
pulse  rectifiers,  the  required  phase  shift  between 
individual  rectifiers  is  given  as  phase  shift  (f>  equals  60° 
divided  by  the  number  of  6-pulse  converters  [16] .  For 
example,  a  12-pulse  converter  would  be  made  from  two  6- 
pulse  converters,  and  60°/2=30°;  therefore,  the  required 
phase  shift  between  the  two  6-pulse  converters  is  30°.  One 
realization  of  such  a  system  is  shown  in  Figure  26.  The 
system  consists  of  two  bridge  rectifiers  connected  so  that 
the  DC  outputs  are  in  series  and  the  two  voltages  add.  The 
AC  inputs  to  the  two  bridges  are  isolated  by  two 
transformers  that  are  connected  wye-delta  and  delta-delta, 
respectively.  The  transformers  are  connected  to  take 
advantage  of  the  differing  phase  shifts  of  the  two 
connections  [ 3 ] . 
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Figure  26.  Series  circuit  arrangement  for  a  12-pulse 

rectifier  using  a  wye-delta  and  delta-delta  connected 
transformer  to  achieve  the  necessary  30°  phase  shift 
between  6-pulse  rectifier  units.  From  [3] . 

While  it  was  relatively  simple  to  realize  the  12-pulse 
rectifier  by  taking  advantage  of  the  30°  phase  shift  that 
naturally  occurs  for  a  wye-delta  transformer,  how  would  we 
go  about  achieving  other  phase  shifts?  For  example,  for  an 
18-  and  24-pulse  count  rectifier,  we  would  require  phase 
shifts  of  20°  and  15°,  respectively.  One  such  device  that 
can  easily  accomplish  this  is  called  the  zig-zag 


transformer  [17]  . 


VI .  SIMULINK®  SIMULATIONS 


In  order  to  validate  the  relationship  between  pulse 
count  and  harmonic  content  (including  THD) ,  simulations  for 
a  6-,  12-,  18-,  and  24-pulse  controlled  rectifier  with  a 

two-pole  LC  output  filter  were  conducted  using  Simulink®, 
specifically  looking  at  the  harmonic  content  of  the  phase  A 
input  line  current  iap  .  Additionally,  it  was  desired  to 

observe  the  behavior  of  the  harmonics  and  %THDj  as  a  was 
varied . 

A.  CIRCUIT  PARAMETERS  COMMON  TO  ALL  SIMULATION  MODELS 

The  fixed  step-size  for  all  simulations  was  10  |a.s 
using  the  Runge-Kutta  (ode4)  solver;  total  simulation  time 
was  1.0  s . 

The  voltage  source  was  a  three-phase  60  Hz  ramp  source 
with  a  2,  980  V/s  slope,  going  from  zero  to  effectively  150 
V  (RMS,  line-to-line)  in  0.05  s.  The  use  of  a  ramp  source 
rather  than  a  step  source  allows  for  faster  steady-state 
results.  This  voltage  source  was  chosen  to  match  that  of  an 
example  problem  from  [4] . 

The  transformers  were  zigzag  phase  shifting 
transformers  with  the  following  common  parameters: 

•  Secondary  winding  (abc)  connection:  Y, 

•  Nominal  power:  1  MVA, 

•  Nominal  frequency:  60  Hz, 

•  Primary  nominal  voltage  Vp  :  150  V  (RMS,  line-to- 

line)  , 

•  Winding  1  resistance  i?,  :  0, 
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•  Winding  1  inductance  Z,  :  0, 

•  Winding  2  resistance  R2  :  0, 

•  Winding  2  inductance  L ,  :  0, 

•  Winding  3  resistance  R3  :  0, 

•  Winding  3  inductance  Z3 :  0, 

•  Magnetizing  branch  resistance  Rm:  10  kf2 ,  and 

•  Magnetizing  branch  inductance  Lm  :  "inf." 

The  main  circuit  component  values  were  as  follows: 

•  Inductor:  23.375  mH, 

•  Capacitor:  200  |iF ,  and 

•  Resistor:  19.3  Q. 

These  values  were  also  chosen  to  match  those  of  the  same 
example  problem  from  [4] . 

For  the  discrete  synchronized  6-pulse  generators,  a 
pulse  width  of  10°  was  used  with  double-pulsing  disabled. 

For  the  6-pulse  controlled  rectifier  module,  an 
inductance  of  1  pH  was  added  to  each  branch  feeding  into 
the  universal  bridges. 

The  universal  bridges  had  the  following  common 
parameters : 

•  Number  of  bridge  arms:  3, 

•  Snubber  resistance  :  100  Q, 

•  Snubber  capacitance  Cs  :  0.1  |iF , 

•  Power  electronic  device:  thyristors, 

•  On-resistance  Ron  :  1  mQ  , 

•  On-inductance  Lon:  0,  and 

•  Forward  voltage  Vf  :  1  mV. 
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The  Simulink®  schematic  for  the  internal  view  of  the 
6-pulse  controlled  rectifier  module  common  to  all  of  the  p- 
pulse  simulation  models  is  shown  in  Figure  27. 
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Figure  27.  Simulink®  schematic  of  internal  view  of 
6-pulse  controlled  rectifier  module. 

B.  DATA  COLLECTION 

As  stated  previously,  the  data  of  interest  was  the 
harmonic  content  of  i  as  a  varied.  Because  it  was 
desirable  to  observe  how  individual  harmonics  (normalized 
to  the  fundamental)  and  %THDj  varied  with  respect  to  a,  it 
was  decided  to  evaluate  the  first  ten  non-zero  harmonics 
(not  including  the  fundamental)  for  each  p-pulse  rectifier. 
The  first  ten  non-zero  harmonics  for  the  simulated  p-pulse 
rectifiers  are  shown  in  Table  9. 
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Table  9.  Harmonic  numbers  for  the  first  ten  non-zero 

harmonics  (not  including  the  fundamental)  for  some 

p-pulse  rectifiers. 


6-pulse 

12-pulse 

18-pulse 

24-pulse 

5 

11 

17 

23 

7 

13 

19 

25 

11 

23 

35 

47 

13 

25 

37 

49 

17 

35 

53 

71 

19 

37 

55 

73 

23 

47 

71 

95 

25 

49 

73 

97 

29 

59 

89 

119 

31 

61 

91 

121 

In  order  to  obtain  the  RMS  values  of  the  individual 
harmonics,  a  string  of  eleven  Fourier  analyser  blocks  was 

used  in  series  with  gain  blocks  set  at  V2/2,  specifying  the 
fundamental  frequency  (60  Hz)  and  harmonic  number  of 
interest  for  each  Fourier  analyser  block.  Additionally,  to 
obtain  %THDi ,  a  total  harmonic  distortion  block  in  series 
with  a  gain  block  equal  to  100  was  used,  again  specifying 
the  fundamental  frequency  of  60  Hz.  The  Simulink®  schematic 
symbology  depicting  this  is  shown  in  Figure  28. 
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Figure  28.  Obtaining  the  RMS  values  for  the 
fundamental  and  first  ten  non-zero  harmonics  and  %THDj 

from  Simulink®. 


Because  it  was  initially  thought  that 
and  %THDj  would  change  faster  at  lower 
(erroneously,  in  this  particular  case) ,  the 
size  for  a  was  chosen  to  be  3°,  starting  at 


the  harmonics 
values  of  a 
initial  Step¬ 
an  0°  and  going 
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up  to  a- 30°,  then  increasing  step-size  to  10°,  and 
continuing  up  to  a  =  80°.  It  was  not  possible  to  go  up  to 
a- 90°  (or  even  a  =  80°  for  all  controlled  rectifiers)  since 
Lcrit  would  approach  infinity,  and  the  controlled  rectifier 
would  no  longer  be  in  CCM. 

The  ten  harmonics  for  each  a  were  then  normalized  to 
the  fundamental  by  dividing  the  RMS  value  of  each  harmonic 
by  the  RMS  value  of  the  fundamental.  The  normalized 
harmonics  were  then  plotted  versus  a  for  each  p-pulse 
controlled  rectifier.  Additionally,  the  %THD:  for  all  four 
of  the  p-pulse  controlled  rectifiers  was  plotted  versus  a  . 

C.  CONTROLLED  6-PULSE  RECTIFIER  SIMULATION 

The  controlled  6-pulse  rectifier  Simulink®  model  is 
shown  in  Figure  29. 


Figure  29.  Simulink®  model  for  a  controlled  6-pulse 
rectifier  with  a  two-pole  LC  output  filter. 

Although  a  transformer  is  not  necessary  for  the  controlled 
6-pulse  rectifier  (because  no  other  phase  shifts  beyond 
those  that  already  exist  for  the  three-phase  source  are 
required) ,  a  zigzag  phase  shifting  transformer  was  used  to 
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provide  isolation.  The  secondary  nominal  voltage  V3  and 
phase  shift  (j)  were  150  V  (RMS,  line-to-line)  and  1X10-6 
degrees  (because  0°  is  not  permitted),  respectively. 

The  first  ten  normalized  harmonics  are  plotted  versus 
a  in  Figure  30.  The  6-pulse  controlled  rectifier  did  not 
have  enough  inductance  at  a  =  80°  to  remain  in  CCM. 
Interestingly,  the  first  normalized  non-zero  harmonic  (h5) 

increases  with  a  (which  causes  %THDi  to  worsen)  ,  but  the 
second  non-zero  harmonic  (h7)  decreases  with  a;  the 
remaining  eight  normalized  harmonics  appear  to  only 
decrease  very  slightly  as  a  increases. 


Figure  30.  First  ten  non-zero  normalized  harmonics 

versus  a  for  a  6-pulse  controlled  rectifier  with  a 
two-pole  LC  output  filter. 
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D.  CONTROLLED  12 -PULSE  RECTIFIER  SIMULATION 

The  controlled  12-pulse  rectifier  Simulink®  model  is 
shown  in  Figure  31. 


6-pulse  controlled 
rectifier  2 

Figure  31.  Simulink®  model  for  a  controlled  12-pulse 
rectifier  with  a  two-pole  LC  output  filter. 

The  secondary  nominal  voltage  V3  for  both  the  upper  and 
lower  zigzag  phase  shifting  transformer  were  75  V  (RMS, 
line-to-line)  while  the  phase  shift  $  was  15°  for  the  upper 
transformer  and  -15°  for  the  lower  transformer.  The  reason 
the  secondary  nominal  voltages  are  the  same  and  are  half  of 
what  they  were  for  the  6-pulse  case  is  because  it  is 
desirable  to  maintain  the  same  overall  rectifier  output 
side  voltage.  Recall  that  since  the  two  6-pulse  controlled 
rectifiers  are  connected  in  series,  their  output  voltages 
are  additive.  The  ±15°  phase  shifts  create  the  necessary  30° 
phase  difference  between  the  two  6-pulse  controlled 
rectifiers  to  ensure  proper  harmonic  cancellation. 
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The  first  ten  normalized  harmonics  are  plotted  versus 
a  in  Figure  32.  The  normalized  harmonics  of  the  12-pulse 
controlled  rectifier  displayed  fairly  similar  behavior  to 
those  of  the  6-pulse  controlled  rectifier,  although,  they 
are  smaller,  as  predicted  by  (177)  and  (178) . 


0  10  20  30  40  50  60  70  80 


a  (degrees) 

— ♦—  hll  -B-hl3  -±-h23  -*-h25  -*-h35 
-#-h37  h47  h49  h59— ♦— h61 


Figure  32.  First  ten  non-zero  normalized  harmonics 

versus  a  for  a  12-pulse  controlled  rectifier  with  a 
two-pole  LC  output  filter. 


E.  CONTROLLED  18-PULSE  RECTIFIER  SIMULATION 

The  controlled  18-pulse  rectifier  Simulink®  model  is 
shown  in  Figure  33. 
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6-pulse  controlled 
rectifier  3 


Figure  33.  Simulink®  model  for  a  controlled  18-pulse 
rectifier  with  a  two-pole  LC  output  filter. 

The  secondary  nominal  voltage  V3  for  all  three  zigzag 
phase  shifting  transformer  was  50  V  (RMS,  line-to-line) 
while  the  phase  shift  ^  was  20°  for  the  upper  transformer, 
1X10-6  degrees  (because  0°  is  not  permitted)  for  the  middle 
transformer,  and  -20°  for  the  lower  transformer.  The  reason 
the  secondary  nominal  voltages  are  the  same  is  for  the  same 
reason  as  before.  The  ±20°  and  0°  (effectively)  phase  shifts 
create  the  necessary  20°  phase  differences  between  the 
three  6-pulse  controlled  rectifiers  to  ensure  proper 
harmonic  cancellation. 

The  first  ten  normalized  harmonics  are  plotted  versus 
a  in  Figure  34.  The  normalized  harmonics  of  the  18-pulse 
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controlled  rectifier  displayed  similar  behavior  to  those  of 
the  other  two  controlled  rectifiers,  although,  they  are 
smaller,  again  as  predicted. 


Figure  34.  First  ten  non-zero  normalized  harmonics 

versus  a  for  an  18-pulse  controlled  rectifier  with  a 
two-pole  LC  output  filter. 


F.  CONTROLLED  24-PULSE  RECTIFIER  SIMULATION 

The  controlled  24-pulse  rectifier  Simulink®  model  is 
shown  in  Figure  35. 
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6-pike  coitroiied 


Figure  35.  Simulink®  model  for  a  controlled  24-pulse 
rectifier  with  a  two-pole  LC  output  filter. 

The  secondary  nominal  voltage  V3  for  all  four  zigzag 
phase  shifting  transformer  was  37.5  V  (RMS,  line-to-line) 
while  the  phase  shift  ^  was  22.5°  for  the  top  transformer, 
7.5°  for  the  second  transformer,  -7.5°  for  the  third 
transformer,  and  -22.5°  for  the  bottom  transformer.  The 
phase  shifts  of  +22.5°  and  +7.5°  create  the  necessary  15° 
phase  differences  between  the  four  6-pulse  controlled 
rectifiers  to  ensure  proper  harmonic  cancellation. 
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The  first  ten  normalized  harmonics  are  plotted  versus 
a  in  Figure  36.  It  is  no  surprise  that  the  normalized 
harmonics  of  the  24-pulse  controlled  rectifier  displayed 
similar  behavior  to  those  of  the  other  three  controlled 
rectifiers,  although,  they  are  once  again  smaller,  just  as 
predicted . 
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Figure  36.  First  ten  non-zero  normalized  harmonics 

versus  a  for  a  24-pulse  controlled  rectifier  with  a 
two-pole  LC  output  filter. 

G.  SIMULATION  RESULTS  FOR  %THDi 

The  %THDj  versus  a  for  p-pulse  controlled  rectifiers 
with  a  two-pole  LC  output  filter  is  shown  in  Figure  37. 
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Figure  37.  Percent  total  harmonic  distortion  in  iap 

versus  a  for  p-pulse  controlled  rectifiers  with  a  two- 

pole  LC  output  filter. 


The  data  clearly  show  that  a  higher  pulse  count  produces  a 
lower  %THDj ,  with  the  24-pulse  controlled  rectifier  the 
clear  winner.  However,  an  incorrect  conclusion  that  might 
be  drawn  from  Figure  37  is  that  %THDi  does  not  vary  much 
with  a  ;  in  fact,  it  appears  as  if  %THDj  only  gets  slightly 
worse  for  very  high  values  of  a  .  The  real  reason  that 
%THDj  does  not  appreciably  worsen  as  a  increases  is 
because  the  value  of  the  inductor  chosen  for  the 
simulations  was  very  large  (23.375  mH)  .  The  point  is  that 
%THDj  and  Lcrit  are  both  related  to  a  in  much  the  same  way. 

That  is.  Figure  37  is  related  to  Figure  15.  If  we  had 
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instead  initially  selected  a  much  smaller  inductor  that 
would  still  ensure  CCM  for  the  6-pulse  controlled  rectifier 
up  to  a- 30°  (  =  3  mH)  ,  we  would  more  than  likely  see  %THDi 

increase  for  much  smaller  values  of  a  ,  particularly  for  the 
6-pulse  controlled  rectifier  and  less  so  for  the  other 
higher  pulse-count  controlled  rectifiers  because  of  the 
reduction  in  Lcrit  as  pulse  count  increases. 
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VII .  CONCLUSIONS 


A.  SUMMARY  OF  RESULTS 

1 .  Critical  Inductance 

The  critical  inductance  as  a  function  of  a  was 
calculated  for  a  p-pulse  controlled  rectifier  with  a  two- 
pole  LC  output  filter.  Plotting  Lcrjt  for  several  multi-pulse 
controlled  rectifiers  demonstrated  how  it  increased  as  a 
increased  and  that  higher  pulse  count  controlled  rectifiers 
have  a  lower  Lcrit  than  lower  pulse  count  controlled 
rectifiers  for  a  given  value  of  a  .  Additionally,  Lcrit  as  a 
percentage  of  Lcrit  for  a  6-pulse  controlled  rectifier  was 
calculated  for  selected  values  of  a  (0°  to  90°,  in  10° 
increments)  for  a  6-,  12-,  18-,  24-,  and  96-pulse 

controlled  rectifier.  This  clearly  illustrated  how  the 
greatest  reduction  in  Lcrit  occurs  in  increasing  from  a  6- 
pulse  to  a  12-pulse  controlled  rectifier  and  that,  while 
further  increasing  pulse  count  does  reduce  Lcrit ,  the 
reduction  is  less. 

2 .  DC  Output  Voltage  Ripple 

The  peak-to-peak  ripple  voltage  as  a  percentage  of  the 
average  output  voltage  as  a  function  of  a  was  calculated 
for  a  p-pulse  controlled  rectifier  with  a  purely  resistive 
load.  Plotting  %Vpp/avg  for  several  multi-pulse  controlled 

rectifiers  demonstrated  how  it  increased  as  a  increased  and 
that  higher  pulse  count  controlled  rectifiers  have  a  lower 
o/oVpp/avg  than  lower  pulse  count  controlled  rectifiers  for  a 
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given  value  of  a  .  Additionally,  %Vpp/avg  as  a  percentage  of 
%Vpp/avg  f°r  a  6-pulse  controlled  rectifier  was  calculated  for 

selected  values  of  a  (0°  to  90°  +  ,  in  10°  increments)  for  a 
6-,  12-,  18-,  24-,  and  96-pulse  controlled  rectifier.  This 

clearly  illustrated  how  the  greatest  reduction  in  %Vf 

occurs  in  increasing  from  a  6-pulse  to  a  12-pulse 
controlled  rectifier  and  that,  while  further  increasing 
pulse  count  does  reduce  %Vpp/avg ,  the  reduction  is  less. 

The  peak-to-peak  ripple  voltage  as  a  percentage  of  the 
average  output  voltage  as  a  function  of  a  was  also 
calculated  for  a  p-pulse  controlled  rectifier  with  a  two- 
pole  LC  output  filter.  Plotting  %Vpp/urg  for  several  multi¬ 
pulse  controlled  rectifiers  demonstrated  how  it  increased 
as  a  increased  and  that  higher  pulse  count  controlled 
rectifiers  have  a  lower  °/oVpplavg  than  lower  pulse  count 
controlled  rectifiers  for  a  given  value  of  a  .  Additionally, 
%VpP/avg  as  a  percentage  of  %V  /avg  for  a  6-pulse  controlled 

rectifier  was  calculated  for  selected  values  of  a  (0°  to 
90°,  in  10°  increments)  for  a  6-,  12-,  18-,  24-,  and  96- 

pulse  controlled  rectifier.  This  clearly  illustrated  how 
the  greatest  reduction  in  %Vpp/avg  occurs  in  increasing  from  a 
6-pulse  to  a  12-pulse  controlled  rectifier  and  that,  while 
further  increasing  pulse  count  does  reduce  %Vt  ,  the 
reduction  is  less. 

3.  Harmonics  and  Total  Harmonic  Distortion 

Simulations  were  conducted  using  Simulink®  for  a  6-, 

12-,  18-,  and  24-pulse  controlled  rectifier  with  a  two-pole 

LC  output  filter,  specifically  investigating  the  behavior 
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of  the  first  ten  non-zero  harmonics  (normalized  to  the 
fundamental)  and  %THDI  of  iup  for  specific  values  of  a  . 

Plotting  the  harmonics  of  each  p-pulse  controlled  rectifier 
validated  that  the  harmonics  for  higher  pulse  count 
controlled  rectifiers  are  smaller  than  those  for  lower 
pulse  count  controlled  rectifiers.  Additionally,  it 
revealed  that  the  first  harmonic  increased  as  a  increased, 
while  the  second  harmonic  decreased  as  a  increased.  The 
remaining  eight  harmonics  did  not  show  significant  change 
as  a  varied. 

Plotting  %THDj  for  the  p-pulse  controlled  rectifiers 
confirmed  that  it  increased  as  a  increased  (though  it  was 
reasoned  that  it  would  show  a  more  definitive  increase  as  a 
increased  if  a  smaller  inductor  had  been  used)  and  that 
higher  pulse  count  controlled  rectifiers  have  a  lower  %THDi 
than  lower  pulse  count  controlled  rectifiers  for  a  given 
value  of  a  .  This  also  illustrated  how  the  greatest 
reduction  in  %THDj  occurs  in  increasing  from  a  6-pulse  to  a 
12-pulse  controlled  rectifier  and  that,  while  further 
increasing  pulse  count  does  reduce  %THDj ,  it  results  in 
diminishing  returns. 

B.  RECOMMENDATIONS  FOR  FUTURE  WORK 

1 .  Discontinuous  Conduction  Mode 

When  the  inductor  in  a  multi-pulse  phase  controlled 
rectifier  with  a  two-pole  LC  output  filter  is  less  than  the 
value  of  Lcrit ,  the  inductor  current  is  no  longer  continuous, 
and  the  equations,  tables,  and  figures  of  Chapters  III  and 
IV  no  longer  apply.  Although  operation  of  a  controlled 
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rectifier  in  discontinuous  conduction  mode  is  not  generally 
considered  desirable,  it  is  still  a  possible  mode  for  the 
rectifier,  and  the  model  is  not  complete  without  it. 

2 .  Size  and  Weight  Considerations 

While  it  is  clear  that  a  higher  pulse  count  controlled 
rectifier  produces  better  results  than  one  with  a  lower 
pulse  count,  it  can  only  do  so  at  the  cost  of  greater 
system  complexity  (and  hence,  less  system  reliability)  and 
with  more  power  electronic  components.  Power  SCRs, 
inductors,  capacitors,  and  especially  low  frequency 
transformers  are  bulky  and  heavy;  this  must  be  considered 
in  real  systems  designed  to  operate  on  Navy  ships  with 
space  and  weight  limitations.  One  useful  study  would  be  to 
use  actual  space  and  weight  characteristics  of  the  power 
electronic  components  for  these  controlled  rectifiers  and 
to  determine  the  optimal  pulse  count  system  that  would 
satisfy  the  electrical  requirements  and  be  within  size  and 
weight  constraints  for  the  intended  Navy  ship  platform. 

3.  Parallel  Vice  Series  Connection  of  Multi-Pulse 
Controlled  Rectifiers 

Only  the  series  (vice  parallel)  connection  of  multiple 
controlled  rectifiers  in  order  to  increase  pulse  count  was 
considered  in  this  thesis.  With  a  series  connection,  the 
output  voltage  of  the  system  is  the  summation  of  the 
individual  output  voltages  of  each  rectifier,  whereas,  with 
a  parallel  connection,  the  output  current  of  the  system  is 
the  summation  of  the  individual  output  currents  of  each 
rectifier.  However,  connecting  rectifiers  in  parallel  is 


more  complex  due  to  fluctuations  in  DC  voltages  which  cause 

circulating  currents  to  flow.  To  remedy  this,  these  types 
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of  systems  typically  require  the  use  of  an  interphase 
transformer  which  absorbs  the  differences  in  the  DC 
voltages  between  individual  rectifiers.  A  useful  study 
would  be  the  evaluation  of  differences  between  the  two 
topologies  as  a  result  of  the  added  interphase  transformer. 

4.  The  Effect  of  a  on  Total  Harmonic  Distortion 


In  general,  increasing  a 

will 

cause 

the  THD  of 

a 

controlled 

rectifier  system 

to 

worsen . 

However, 

a 

theoretical 

mathematical  model 

for 

how  THD 

behaves  as 

a 

function  of  a  similar  to  how  Lcrjt  behaves  as  a  changes  for 
a  p-pulse  controlled  rectifier  with  a  two-pole  LC  output 
filter  (as  in  Chapter  III)  would  be  beneficial. 

5 .  A  Better  Solution  for  the  Output  Voltage  Ripple 
for  a  Controlled  Rectifier  with  a  Two-Pole  LC 
Output  Filter 

In  Chapter  IV,  approximations  were  introduced  in  the 
determination  of  the  output  voltage  ripple  for  a  controlled 
rectifier  with  a  two-pole  LC  output  filter.  This  had  to  do 

with  the  fact  that  \lAc\  was  calculated  from  VAC  and  the 

overall  load  impedance  |Ztoto/|  only  at  the  frequency  pco  .  While 
po)  is  the  dominant  frequency,  load  impedances  at  frequency 
multiples  of  pco  also  exist.  A  more  accurate  estimate  of 

\Iac\  can  be  obtained  by  using  the  principle  of  superposition 

to  solve  for  ltd  using  several  values  of  \Ztotal\  at  frequency 
multiples  of  pco  and  adding  them  together.  While  this  will 
provide  a  better  estimate  of  \lAc\  '  which  will  in  turn  result 
in  a  more  accurate  value  for  VC(RMS) ,  the  fact  remains  that 
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Vc(rms)  is  not  perfectly  sinusoidal,  and  VC(pp_ripple)  and  %Vc{pplavg) 

will  still  be  approximations,  although,  somewhat  better 
than  before.  It  would  still  be  interesting  to  determine  the 
error  caused  by  this  approximation. 
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APPENDIX.  MATLAB®  CODE 


A.  BASIC  LINE-TO-LINE  VOLTAGE  WAVEFORM  PLOT  FOR  A 
HYPOTHETICAL  THREE-PHASE  CIRCUIT 

clc 

clear 

f  =  6  0  ; 
w=2*pi*f ; 

t=linspace ( 0 ,  3/pi,  1000); 

Vab  =  sind(w*t); 

Vbc  =  sind (w*t-120) ; 

Vca  =  sind (w*t-240) ; 

Vba=-Vab; 

Vcb=-Vbc; 

Vac=-Vca; 

hold  on 

plot(w*t,  Vab,  'b'  ,  w*t,  Vbc,  'g'  ,  w*t,  Vca,  'r'  ,  'LineWidth' ,  2) 

axis ( [0,360, -1,1] ) 

xlabel ( 'angle  (degrees)') 

ylabel ( 'Voltage  (Volts)') 

legend ( 'V_a_b' ,  'V_b_c' ,  'V_c_a' ) 

grid  on 

hold  off 

B.  LINE-TO-LINE  VOLTAGE  WAVEFORM  PLOT  SHOWING  SIX  PULSES 
FOR  A  HYPOTHETICAL  THREE-PHASE  CIRCUIT 


clc 

clear 

f  =  6  0  ; 
w=2*pi*f  ; 

t=linspace ( 0 ,  3/pi,  1000); 

Vab  =  sind(w*t); 

Vbc  =  sind (w*t-120) ; 

Vca  =  sind (w*t-240) ; 

Vba=-Vab; 

Vcb=-Vbc; 

Vac=-Vca; 

hold  on 

plot(w*t,  Vab,  'b' ,  w*t,  Vbc,  'g' ,  w*t,  Vca,  'r' ,  'LineWidth',  2) 
plot(w*t,  Vba,  'c' ,  w*t,  Vcb,  'k'  ,  w*t.  Vac,  'm' ,  'LineWidth',  2) 
axis ( [0, 360, -1,1]) 
xlabel ( 'angle  (degrees)') 
ylabel ( 'Voltage  (Volts)') 

legend ( 'V_a_b' ,  'V_b_c' ,  'V_c_a' ,  'V_b_a' ,  'V_c_b' ,  'V_a_c' ) 

grid  on 
hold  off 
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C.  CRITICAL  INDUCTANCE  CALCULATION  AND  PLOT  FOR  A 
CONTROLLED  6-PULSE  RECTIFIER 


clear 

clc 

p  =  6;  %  p  =  pulse  count;  i.e.,  p  =  6  for  6-pulse  controlled  rectifier 
lil  =  90- (360/ (2*p) ) ;  %  lil  =  lower  integration  limit  for  start  of  a 
pulse 

A  =  cosd(lil);  %  constant  to  make  math  cleaner/neater 
B  =  sind(lil);  %  constant  to  make  math  cleaner/neater 
alpha  u  =  85;  %  upper  limit  in  degrees  of  firing  angle  alpha 
num  points  =  5e5;  %  number  of  points  for  alpha 

alpha  b  r  =  atan ( (p/pi) - (B/A) ) ;  %  boundary  angle  alpha  in  radians  that 
%  defines  transition  point  between  submodes  of  continuous 
conduction 

alpha  b  d  =  alpha  b  r*180/pi;  %  boundary  angle  alpha  in  degrees 
w  =  2*pi*60;  %  radian  frequency  for  60Hz 
R  =  1;  %  unit  resistance  chosen  for  ease  of  comparison 

%  code  for  angles  greater  than  the  boundary  angle 

alpha  d  =  linspace (alpha  b  d,  alpha  u,  num  points);  %  range  of  alpha 
%  in  degrees  for  Lcritl 

Lcritl  =  ( (R*tand (alpha_d) ) /w) * ( 1- ( (pi*B) / (p*A) ) ) ;  %  critical 
inductance 

%  for  angles  greater  than  the  boundary  angle 
hold  on 

plot (alpha  d,  Lcritl,  'k-.',  'LineWidth' ,  2) 

%  code  for  angles  less  than  the  boundary  angle 

alpha  d  =  linspace (0,  alpha  b  d,  num  points);  %  range  of  alpha 
%  in  degrees  for  Lcrit2 
alpha  r  =  alpha  d*pi/180;  %  alpha  in  radians 

theta_d  =  asind ( (p*A*cosd (alpha_d) ) /pi)  -  lil;  %  angle  in  degrees  at 
which 

%  the  inductor  current  will  be  minimum  for  this  submode 
theta  r  =  theta  d*pi/180;  %  theta  in  radians 
Lcrit2  =  (R. / (p*A*w. *cosd (alpha  d) ) ) . * . . . 

(pi . *cosd (theta_d+lil)  +  . .  . 

(p*A. *sind (alpha_d) ) + . . . 

( (p. *theta_r-p. *alpha_r-pi) *A. *cosd (alpha_d) ) ) ; 
plot (alpha  d,  Lcrit2,  'k' ,  'LineWidth',  2) 

%  title ( 'Critical  Inductance  vs.  \alpha  for  a  6  Pulse  Controlled 
Rectifier'  ) 

xlabel ( ' \alpha  (degrees)') 

ylabel (  'Critical  Inductance  (Henry) '  ) 

legend  (' \alpha  \geq  10 . 08 A\circ' ,  '/alpha  \leq  10 . 08  A\circ' ,  .  .  . 

'location',  'North') 
grid  on 
hold  off 
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D.  CRITICAL  INDUCTANCE  CALCULATIONS  AND  PLOT  FOR 
CONTROLLED  6-,  12- ,  18- ,  24-,  AND  96-PULSE  RECTIFIERS 


clear 

clc 

hold  on 

alpha  u  =  30;  %  upper  limit  in  degrees  of  firing  angle  alpha 
num  points  =  5e5;  %  number  of  points  for  alpha 

%%  6  pulse 

p  =  6;  %  p  =  pulse  count;  i.e.,  p  =  6  for  6-pulse  controlled  rectifier 
lil  =  90- (360/ (2*p) ) ;  %  lil  =  lower  integration  limit  for  start  of  a 
pulse 

A  =  cosd(lil);  %  constant  to  make  math  cleaner/neater 
B  =  sind(lil);  %  constant  to  make  math  cleaner/neater 

alpha  b  r  =  atan ( (p/pi) - (B/A) ) ;  %  boundary  angle  alpha  in  radians  that 
%  defines  transition  point  between  submodes  of  continuous 
conduction 

alpha  b  d  =  alpha  b  r*180/pi;  %  boundary  angle  alpha  in  degrees 
w  =  2*pi*60;  %  radian  frequency  for  60Hz 
R  =  1;  %  unit  resistance  chosen  for  ease  of  comparison 

%  code  for  angles  greater  than  the  boundary  angle 

alpha  d  =  linspace (alpha  b  d,  alpha  u,  num  points);  %  range  of  alpha 
%  in  degrees  for  Lcritl 

Lcritl  =  ( (R*tand (alpha_d) ) /w) * (1- ( (pi*B) / (p*A) ) ) ;  %  critical 
inductance 

%  for  angles  greater  than  the  boundary  angle 
plot (alpha  d,  Lcritl,  'r-.',  'LineWidth' ,  2) 

%  code  for  angles  less  than  the  boundary  angle 

alpha  d  =  linspace (0,  alpha  b  d,  num  points);  %  range  of  alpha 
%  in  degrees  for  Lcrit2 
alpha  r  =  alpha  d*pi/180;  %  alpha  in  radians 

theta_d  =  asind ( (p*A*cosd (alpha_d) ) /pi)  -  lil;  %  angle  in  degrees  at 
which 

%  the  inductor  current  will  be  minimum  for  this  submode 
theta_r  =  theta_d*pi/180;  %  theta  in  radians 
Lcrit2  =  (R. / (p*A*w. *cosd (alpha_d) )).*.. . 

(pi . *cosd (theta_d+lil)  +  . .  . 

(p*A. *sind (alpha_d) ) + . . . 

( (p. *theta_r-p. *alpha_r-pi) *A. *cosd (alpha_d) ) ) ; 
plot (alpha  d,  Lcrit2,  'r' ,  'LineWidth',  2) 

%%  12  pulse 

p  =  12;  %  p  =  pulse  count;  i.e.,  p  =  6  for  6-pulse  controlled  rectifier 
lil  =  90- (360/ (2*p) ) ;  %  lil  =  lower  integration  limit  for  start  of  a 
pulse 

A  =  cosd(lil);  %  constant  to  make  math  cleaner/neater 
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B  =  sind(lil);  %  constant  to  make  math  cleaner/neater 

alpha  b  r  =  atan ( (p/pi) - (B/A) ) ;  %  boundary  angle  alpha  in  radians  that 
%  defines  transition  point  between  submodes  of  continuous 
conduction 

alpha  b  d  =  alpha  b  r*180/pi;  %  boundary  angle  alpha  in  degrees 
w  =  2*pi*60;  %  radian  frequency  for  60Hz 
R  =  1;  %  unit  resistance  chosen  for  ease  of  comparison 

%  code  for  angles  greater  than  the  boundary  angle 

alpha  d  =  linspace (alpha  b  d,  alpha  u,  num  points);  %  range  of  alpha 
%  in  degrees  for  Lcritl 

Lcrit3  =  ( (R*tand (alpha_d) ) /w) * (1- ( (pi*B) / (p*A) ) ) ;  %  critical 
inductance 

%  for  angles  greater  than  the  boundary  angle 
plot (alpha  d,  Lcrit3,  'm- . ' ,  'LineWidth' ,  2) 

%  code  for  angles  less  than  the  boundary  angle 

alpha  d  =  linspace (0,  alpha  b  d,  num  points);  %  range  of  alpha 
%  in  degrees  for  Lcrit2 
alpha  r  =  alpha  d*pi/180;  %  alpha  in  radians 

theta_d  =  asind ( (p*A*cosd (alpha_d) ) /pi)  -  lil;  %  angle  in  degrees  at 
which 

%  the  inductor  current  will  be  minimum  for  this  submode 
theta  r  =  theta  d*pi/180;  %  theta  in  radians 
Lcrit4  =  (R. / (p*A*w. *cosd (alpha  d) ) ) . * . . . 

(pi . *cosd (theta_d+lil)  +  . .  . 

(p*A. *sind (alpha_d) ) + . . . 

( (p. *theta_r-p. *alpha_r-pi) *A. *cosd (alpha_d) ) ) ; 
plot (alpha  d,  Lcrit4,  'm' ,  'LineWidth',  2) 

%%  18  pulse 

p  =  18;  %  p  =  pulse  count;  i.e.,  p  =  6  for  6-pulse  controlled  rectifier 
lil  =  90- (360/ (2*p) ) ;  %  lil  =  lower  integration  limit  for  start  of  a 
pulse 

A  =  cosd(lil);  %  constant  to  make  math  cleaner/neater 
B  =  sind(lil);  %  constant  to  make  math  cleaner/neater 

alpha  b  r  =  atan ( (p/pi) - (B/A) ) ;  %  boundary  angle  alpha  in  radians  that 
%  defines  transition  point  between  submodes  of  continuous 
conduction 

alpha  b  d  =  alpha  b  r*180/pi;  %  boundary  angle  alpha  in  degrees 
w  =  2*pi*60;  %  radian  frequency  for  60Hz 
R  =  1;  %  unit  resistance  chosen  for  ease  of  comparison 

%  code  for  angles  greater  than  the  boundary  angle 

alpha  d  =  linspace (alpha  b  d,  alpha  u,  num  points);  %  range  of  alpha 
%  in  degrees  for  Lcritl 

Lcrit5  =  ( (R*tand (alpha_d) ) /w) * ( 1- ( (pi*B) / (p*A) ) ) ;  %  critical 
inductance 

%  for  angles  greater  than  the  boundary  angle 
plot (alpha  d,  Lcrit5,  'g- . ' ,  'LineWidth',  2) 

%  code  for  angles  less  than  the  boundary  angle 

alpha  d  =  linspace (0,  alpha  b  d,  num  points);  %  range  of  alpha 
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angle  in  degrees  at 


%  in  degrees  for  Lcrit2 
alpha  r  =  alpha  d*pi/180;  %  alpha  in  radians 
theta_d  =  asind ( (p*A*cosd (alpha_d) ) /pi)  -  lil;  % 
which 

%  the  inductor  current  will  be  minimum  for  this  submode 
theta_r  =  theta_d*pi/180;  %  theta  in  radians 
Lcrit6  =  (R. / (p*A*w. *cosd (alpha_d) )).*.. . 

(pi . *cosd (theta_d+lil) + . . . 

(p*A. *sind (alpha_d) ) + . . . 

( (p. *theta_r-p. *alpha_r-pi) *A. *cosd (alpha_d) ) ) ; 
plot (alpha  d,  Lcrit6,  'g' ,  'LineWidth' ,  2) 

%%  24  pulse 

p  =  24;  %  p  =  pulse  count;  i.e.,  p  =  6  for  6-pulse  controlled  rectifier 
lil  =  90- (360/ (2*p) ) ;  %  lil  =  lower  integration  limit  for  start  of  a 
pulse 

A  =  cosd(lil);  %  constant  to  make  math  cleaner/neater 
B  =  sind(lil);  %  constant  to  make  math  cleaner/neater 

alpha  b  r  =  atan ( (p/pi ) - (B/A) ) ;  %  boundary  angle  alpha  in  radians  that 
%  defines  transition  point  between  submodes  of  continuous 
conduction 

alpha  b  d  =  alpha  b  r*180/pi;  %  boundary  angle  alpha  in  degrees 
w  =  2*pi*60;  %  radian  frequency  for  60Hz 
R  =  1;  %  unit  resistance  chosen  for  ease  of  comparison 

%  code  for  angles  greater  than  the  boundary  angle 

alpha  d  =  linspace (alpha  b  d,  alpha  u,  num  points);  %  range  of  alpha 
%  in  degrees  for  Lcritl 

Lcrit7  =  ( (R*tand (alpha_d) ) /w) * (1- ( (pi*B) / (p*A) ) ) ;  %  critical 
inductance 

%  for  angles  greater  than  the  boundary  angle 
plot (alpha  d,  Lcrit7,  'k-.',  'LineWidth',  2) 

%  code  for  angles  less  than  the  boundary  angle 

alpha  d  =  linspace (0,  alpha  b  d,  num  points);  %  range  of  alpha 
%  in  degrees  for  Lcrit2 
alpha  r  =  alpha  d*pi/180;  %  alpha  in  radians 

theta_d  =  asind ( (p*A*cosd (alpha_d) ) /pi)  -  lil;  %  angle  in  degrees  at 
which 

%  the  inductor  current  will  be  minimum  for  this  submode 
theta  r  =  theta  d*pi/180;  %  theta  in  radians 
Lcrit8  =  (R. / (p*A*w. *cosd (alpha_d) )).*.. . 

(pi . *cosd (theta_d+lil) + . . . 

(p*A. *sind (alpha_d) ) + . . . 

( (p. *theta_r-p. *alpha_r-pi) *A. *cosd (alpha_d) ) ) ; 
plot (alpha  d,  Lcrit8,  'k' ,  'LineWidth',  2) 

%%  96  pulse 

p  =  96;  %  p  =  pulse  count;  i.e.,  p  =  6  for  6-pulse  controlled  rectifier 
lil  =  90- (360/ (2*p) ) ;  %  lil  =  lower  integration  limit  for  start  of  a 
pulse 

A  =  cosd(lil);  %  constant  to  make  math  cleaner/neater 
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B  =  sind(lil);  %  constant  to  make  math  cleaner/neater 

alpha  b  r  =  atan ( (p/pi) - (B/A) ) ;  %  boundary  angle  alpha  in  radians  that 
%  defines  transition  point  between  submodes  of  continuous 
conduction 

alpha  b  d  =  alpha  b  r*180/pi;  %  boundary  angle  alpha  in  degrees 
w  =  2*pi*60;  %  radian  frequency  for  60Hz 
R  =  1;  %  unit  resistance  chosen  for  ease  of  comparison 

%  code  for  angles  greater  than  the  boundary  angle 

alpha  d  =  linspace (alpha  b  d,  alpha  u,  num  points);  %  range  of  alpha 
%  in  degrees  for  Lcritl 

Lcrit9  =  ( (R*tand (alpha_d) ) /w) * (1- ( (pi*B) / (p*A) ) ) ;  %  critical 
inductance 

%  for  angles  greater  than  the  boundary  angle 
plot (alpha  d,  Lcrit9,  'b- . ' ,  'LineWidth' ,  2) 

%  code  for  angles  less  than  the  boundary  angle 

alpha  d  =  linspace (0,  alpha  b  d,  num  points);  %  range  of  alpha 
%  in  degrees  for  Lcrit2 
alpha  r  =  alpha  d*pi/180;  %  alpha  in  radians 

theta_d  =  asind ( (p*A*cosd (alpha_d) ) /pi)  -  lil;  %  angle  in  degrees  at 
which 

%  the  inductor  current  will  be  minimum  for  this  submode 
theta  r  =  theta  d*pi/180;  %  theta  in  radians 
LcritlO  =  (R. / (p*A*w. *cosd (alpha  d) ) )  . * .  .  . 

(pi . *cosd (theta_d+lil) + . . . 

(p*A. *sind (alpha_d) ) + . . . 

( (p. *theta_r-p. *alpha_r-pi) *A. *cosd (alpha_d) ) ) ; 
plot (alpha  d,  LcritlO,  'b' ,  'LineWidth',  2) 

%%  Plot  Formatting 

grid  on 

%%title ( 'Critical  Inductance  vs.  \alpha  for  an  n-Pulse  Controlled 
Rectifier' ) 

xlabel ( ' \alpha  (degrees)') 

ylabel ( 'Critical  Inductance  (Henry) '  ) 

legend (' 6-pulse,  \alpha  \geq  10 . 08 A\circ' ,  '6-pulse,  \alpha  \leq 
10 . 08 A\circ'  ,  .  .  . 


\ — 1 
o 

LO 

'12-pulse, 
'\circ' ,  .  .  . 

\alpha 

\geq 

5 . 01 A\circ' , 

'12-pulse,  \alpha  \leq 

3.34' 

'18-pulse, 
'\circ' ,  .  .  . 

\alpha 

\geq 

3 . 34 A\circ' , 

'18-pulse,  \alpha  \leq 

2.50' 

'24-pulse, 
'\circ' ,  .  .  . 

\alpha 

\geq 

2 . 50 A\circ' , 

'24-pulse,  \alpha  \leq 

'96-pulse, 

\alpha 

\geq 

0 . 63A\circ' , 

'96-pulse,  \alpha  \leq 

0 . 63A\circ' , . . . 

'location',  'Northwest') 
hold  off 

E.  VOLTAGE  RIPPLE  CALCULATIONS  AND  PLOT  FOR  A  CONTROLLED 
6-PULSE  RECTIFIER  WITH  A  RESISTIVE  LOAD 


clear 
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clc 


%%  6-pulse 


p  =  6;  %  p  =  pulse  count;  i.e.,  p  =  6  for  6-pulse  controlled  rectifier 
Vm  =  1; 


%  lil  =  lower  integration  limit  for  start  of  a  pulse 
%  uil  =  upper  integration  limit  for  start  of  a  pulse 
constant  to  make  math  cleaner/neater 


lil  =  90- (180/p) ; 
uil  =  90+ (180/p) ; 

A  =  cosd  (lil) ;  % 

E  =  cosd (2*lil) ; 

F  =  sind (2*lil) ; 

alpha  bl  d  =  (uil-lil) /2;  %  first  boundary  angle  in  degrees 
alpha  b2  d  =  lil;  %  second  boundary  angle  in  degrees 
num  points  =  1000;  %  number  of  points  for  alpha 


%  code  for  alpha  <  boundary  angle  1 

alpha  d  =  linspace(0,  alpha  bl  d,  num  points);  %  range  of  alpha  in 
degrees 

Vavgl  =  (A*Vm*p*cosd (alpha  d))/pi; 

Vrmsl  =  (Vm/sqrt (2) ) *sqrt (1  +  ( (F*p*cosd (2*alpha_d) ) / (2*pi) ) ) ; 

Vrms  ripplel  =  sqrt (Vrmsl . A2  -  Vavgl. A2); 

Vpp  ripplel  =  Vm* ( 1-sind (alpha  d  +  uil)); 

Vpp  by  Vavg  percentl  =  (Vpp  ripplel  .  /  Vavgl) *100; 
figure ( 1 ) 

plot (alpha  d,  Vavgl,  'k' ,  .  .  . 
alpha  d,  Vrmsl,  'r' , . . . 
alpha  d,  Vrms  ripplel,  'b' , . . . 
alpha  d,  Vpp  ripplel,  'm' ,  'LineWidth' ,  2) 
xlabel ( ' \alpha  (degrees)') 
ylabel ( 'Voltage  (Volts)') 

legend  (  'Vavg',  'VRMS',  'V  R  M  S  -  r  i  p  p  1  e'  ,  'Vpp- 
ripple',  'location',  'NorthEast' ) 
grid  on 
hold  on 
figure  (2 ) 

plot (alpha  d,  Vpp  by_Vavg  percentl,  'k' ,  'LineWidth',  2) 
xlabel (' \alpha  (degrees)') 
ylabel ( '%V_p_p_/_a_v_g  (%)  '  ) 
grid  on 

axis([0  uil  0  1000]) 
hold  on 


%  code  for  boundary  angle  1  <  alpha  <  boundary  angle  2 

alpha  d  =  linspace (alpha  bl  d,  alpha  b2  d,  num  points);  %  range  of 

alpha  in  degrees 

Vavg2  =  (A*Vm*p*cosd (alpha  d))/pi; 

Vrms2  =  (Vm/sqrt (2 )) *sqrt ( 1  +  ( (F*p*cosd (2*alpha_d) ) / (2*pi) ) ) ; 

Vrms  ripple2  =  sqrt (Vrms2 . A2  -  Vavg2.A2); 

Vpp  ripple2  =  Vm* (sind (alpha  d  +  lil)  -  sind(alpha  d  +  uil)); 

Vpp  by  Vavg  percent2  =  (Vpp  ripple2  .  /  Vavg2)*100; 
figure ( 1 ) 

plot (alpha  d,  Vavg2,  'k' , . . . 
alpha  d,  Vrms2,  'r' ,  .  .  . 
alpha  d,  Vrms  ripple2,  'b' ,  .  .  . 
alpha  d,  Vpp  ripple2,  'm' ,  'LineWidth',  2) 
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figure  (2) 

plot (alpha  d,  Vpp  by  Vavg  percent2,  'k' ,  'LineWidth' ,  2) 


range  of  alpha  in 


%  code  for  alpha  >  boundary  angle  2 
alpha  d  =  linspace (alpha  b2  d,  uil,  num  points); 
degrees 

Vavg3  =  ( (Vm*p) / (2*pi) ) * (1  +  cosd(alpha_d  +  lil)  )  ; 

Vrms3  =  (Vm/sqrt (2 )) *sqrt (( (2+p) /4 )-( (alpha_d*p) /360 )+.. . 

( (E*p*sind (2*alpha_d) ) / (4*pi) )+( (F*p*cosd (2*alpha_d) ) / (4*pi) ) ) ; 
Vrms  ripple3  =  sqrt (Vrms3 . A2  -  Vavg3.A2); 

Vpp  ripple3  =  Vm*sind (alpha  d  +  lil); 

Vpp_by_Vavg_percent3  =  (Vpp_ripple3  ./  Vavg3) *100; 
figure ( 1 ) 

Vavg3,  'k' ,  .  .  . 

Vrms 3,  'r' ,  .  .  . 

Vrms  ripple3,  'b'  ,  .  .  . 

Vpp  ripple3,  'm'  ,  'LineWidth',  2) 


plot (alpha_d, 
alpha_d, 
alpha_d, 
alpha_d, 
figure  (2) 
plot (alpha  d. 


Vpp  by  Vavg  percent3,  'k' ,  'LineWidth',  2) 


hold  off 


F.  VOLTAGE  RIPPLE  CALCULATIONS  AND  PLOT  FOR  A  CONTROLLED 
6-,  12- ,  18- ,  24-,  AND  96-PULSE  RECTIFIER  WITH  A 

RESISTIVE  LOAD 

clear 

clc 


hold  on 
%%  6-pulse 


p  =  pulse  count;  i.e.,  p  =  6  for  6-pulse  controlled  rectifier 


p  =  6; 

Vm  =  1  ; 

lil  =  90- (180/p) 
uil  =  90+ (180/p) 

A  =  cosd (lil) ;  \ 

E  =  cosd (2*lil) ; 

F  =  sind (2*lil) ; 
alpha  bl  d  =  (uil-lil)/2; 
alpha  b2  d  =  lil; 
num  points  =  1000; 


%  lil  =  lower  integration  limit  for  start  of 
%  uil  =  upper  integration  limit  for  start  of 
constant  to  make  math  cleaner/neater 


first  boundary  angle  in  degrees 
second  boundary  angle  in  degrees 
number  of  points  for  alpha 


pulse 

pulse 


%  code  for  alpha  <  boundary  angle  1 

alpha  d  =  linspace (0,  alpha  bl  d,  num  points);  %  range  of  alpha  in 
degrees 

Vavgl  =  (A*Vm*p*cosd (alpha  d))/pi; 

Vrmsl  =  (Vm/sqrt (2 )) *sqrt ( 1  +  ( (F*p*cosd (2*alpha_d) ) / (2*pi) ) ) ; 

Vrms  ripplel  =  sqrt (Vrmsl . A2  -  Vavgl. A2); 

Vpp  ripplel  =  Vm* ( 1-sind (alpha  d  +  uil)); 

Vpp  by  Vavg  percentl  =  (Vpp  ripplel  ./  Vavgl) *100; 
plot (alpha  d,  Vpp  by  Vavg  percentl,  'r' ,  'LineWidth',  2) 


%  %  code  for  boundary  angle  1  <  alpha  <  boundary  angle  2 
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range  of 


%  alpha  d  =  linspace  (alpha  bl_d,  alpha  b2  d,  num_points)  ;  % 
alpha  in  degrees 

%  Vavg2  =  (A*Vm*p*cosd (alpha  d) ) /pi; 

%  Vrms2  =  (Vm/sqrt (2 )) *sqrt ( 1  +  ( (F*p*cosd (2*alpha_d) ) / (2*pi) ) ) ; 

%  Vrms  ripple2  =  sqrt (Vrms2 . A2  -  Vavg2.A2); 

%  Vpp  ripple2  =  Vm* (sind (alpha  d  +  lil)  -  sind(alpha  d  +  uil)); 

%  Vpp  by_Vavg  percent2  =  (Vpp  ripple2  ./  Vavg2)*100; 

%  plot (alpha  d,  Vpp  by_Vavg  percent2,  'r--'  ,  'LineWidth' ,  2) 

O, 

o 

%  %  code  for  alpha  >  boundary  angle  2 

%  alpha  d  =  linspace (alpha  b2_d,  uil,  num  points);  %  range  of  alpha  in 
degrees 

%  Vavg3  =  ( (Vm*p) / (2*pi) ) * (1  +  cosd(alpha_d  +  lil)); 

%  Vrms3  =  (Vm/sqrt (2 )) *sqrt (( (2+p) /4 )-( (alpha_d*p) /360 )+.. . 

O, 

o 

( (E*p*sind (2*alpha_d) ) / (4*pi) ) + ( (F*p*cosd (2*alpha_d) ) / (4*pi) ) ) ; 

%  Vrms  ripple3  =  sqrt (Vrms3 . A2  -  Vavg3.A2); 

%  Vpp  ripple3  =  Vm*sind (alpha_d  +  lil); 

%  Vpp_by_Vavg_percent3  =  (Vpp_ripple3  ./  Vavg3)*100; 

%  plot (alpha  d,  Vpp  by_JVavg  percent3,  'r-.',  'LineWidth',  2) 

%%  12-pulse 

p  =  12;  %  p  =  pulse  count;  i.e.,  p  =  6  for  6-pulse  controlled  rectifier 
Vm  =  1  ; 

lil  =  90- (180/p);  %  lil  =  lower  integration  limit  for  start  of  a  pulse 
uil  =  90+ (180/p);  %  uil  =  upper  integration  limit  for  start  of  a  pulse 
A  =  cosd(lil);  %  constant  to  make  math  cleaner/neater 
E  =  cosd (2*lil) ; 

F  =  sind (2*lil) ; 

alpha  bl  d  =  (uil-lil)/2;  %  first  boundary  angle  in  degrees 
alpha  b2  d  =  lil;  %  second  boundary  angle  in  degrees 
num  points  =  1000;  %  number  of  points  for  alpha 

%  code  for  alpha  <  boundary  angle  1 

alpha  d  =  linspace (0,  alpha  bl  d,  num  points);  %  range  of  alpha  in 
degrees 

Vavgl  =  (A*Vm*p*cosd (alpha  d) ) /pi; 

Vrmsl  =  (Vm/sqrt (2 )) *sqrt ( 1  +  ( (F*p*cosd (2*alpha_d) ) / (2*pi) ) ) ; 

Vrms  ripplel  =  sqrt (Vrmsl . A2  -  Vavgl. A2 ); 

Vpp  ripplel  =  Vm* ( 1-sind (alpha  d  +  uil)); 

Vpp  by  Vavg  percentl  =  (Vpp  ripplel  ./  Vavgl) *100; 
plot (alpha  d,  Vpp  by  Vavg  percentl,  'm' ,  'LineWidth',  2) 

%  code  for  boundary  angle  1  <  alpha  <  boundary  angle  2 

alpha  d  =  linspace (alpha  bl  d,  alpha  b2  d,  num  points);  %  range  of 

alpha  in  degrees 

Vavg2  =  (A*Vm*p*cosd (alpha  d))/pi; 

Vrms2  =  (Vm/sqrt (2 )) *sqrt ( 1  +  ( (F*p*cosd (2*alpha_d) ) / (2*pi) ) ) ; 

Vrms  ripple2  =  sqrt (Vrms2 . A2  -  Vavg2.A2); 

Vpp  ripple2  =  Vm* (sind (alpha  d  +  lil)  -  sind(alpha  d  +  uil)); 

Vpp  by  Vavg  percent2  =  (Vpp  ripple2  ./  Vavg2)*100; 

plot (alpha  d,  Vpp  by_Vavg  percent2,  'm--' ,  'LineWidth',  2) 


code  for  alpha  >  boundary  angle  2 
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range  of  alpha  in 


%  alpha  d  =  linspace  (alpha  b2_d,  uil,  num  points);  % 
degrees 

%  Vavg3  =  ( (Vm*p) / (2*pi) ) * (1  +  cosd(alpha_d  +  lil) ) ; 

%  Vrms3  =  (Vm/sqrt (2 ) ) *sqrt ( ( (2+p) /4 ) - ( (alpha_d*p) /360 ) + . . . 

O, 

o 

( (E*p*sind (2*alpha_d) ) / (4*pi) ) + ( (F*p*cosd (2*alpha_d) ) / (4*pi) ) ) ; 

%  Vrms  ripple3  =  sqrt (Vrms3 . A2  -  Vavg3.A2); 

%  Vpp  ripple3  =  Vm*sind (alpha^d  +  lil)  ; 

%  Vpp_by_Vavg_percent3  =  (Vpp_ripple3  ./  Vavg3)*100; 

%  plot (alpha  d,  Vpp  by_Vavg  percent3,  'm- . ' ,  'LineWidth' ,  2) 

%%  18-pulse 

p  =  18;  %  p  =  pulse  count;  i.e.,  p  =  6  for  6-pulse  controlled  rectifier 
Vm  =  1  ; 

lil  =  90- (180/p);  %  lil  =  lower  integration  limit  for  start  of  a  pulse 
uil  =  90+ (180/p);  %  uil  =  upper  integration  limit  for  start  of  a  pulse 
A  =  cosd(lil);  %  constant  to  make  math  cleaner/neater 
E  =  cosd (2*lil) ; 

F  =  sind (2*lil) ; 

alpha  bl  d  =  (uil-lil)/2;  %  first  boundary  angle  in  degrees 
alpha  b2  d  =  lil;  %  second  boundary  angle  in  degrees 
num  points  =  1000;  %  number  of  points  for  alpha 

%  code  for  alpha  <  boundary  angle  1 

alpha  d  =  linspace (0,  alpha  bl  d,  num  points);  %  range  of  alpha  in 
degrees 

Vavgl  =  (A*Vm*p*cosd (alpha  d) ) /pi; 

Vrmsl  =  (Vm/sqrt (2 )) *sqrt ( 1  +  ( (F*p*cosd (2*alpha_d) ) / (2*pi) ) ) ; 

Vrms  ripplel  =  sqrt (Vrmsl . A2  -  Vavgl. A2); 

Vpp  ripplel  =  Vm* ( 1-sind (alpha  d  +  uil)); 

Vpp  by  Vavg  percentl  =  (Vpp  ripplel  .  /  Vavgl) *100; 
plot (alpha  d,  Vpp  by  Vavg  percentl,  'g' ,  'LineWidth',  2) 

%  code  for  boundary  angle  1  <  alpha  <  boundary  angle  2 

alpha  d  =  linspace (alpha  bl  d,  alpha  b2  d,  num  points);  %  range  of 

alpha  in  degrees 

Vavg2  =  (A*Vm*p*cosd (alpha  d))/pi; 

Vrms2  =  (Vm/sqrt (2 )) *sqrt ( 1  +  ( (F*p*cosd (2*alpha_d) ) / (2*pi) ) ) ; 

Vrms  ripple2  =  sqrt (Vrms2 . A2  -  Vavg2.A2); 

Vpp  ripple2  =  Vm* (sind (alpha  d  +  lil)  -  sind(alpha  d  +  uil)); 

Vpp  by  Vavg  percent2  =  (Vpp  ripple2  ./  Vavg2)*100; 

plot (alpha  d,  Vpp  by_Vavg  percent2,  'g--' ,  'LineWidth',  2) 

%  %  code  for  alpha  >  boundary  angle  2 

%  alpha  d  =  linspace (alpha  b2  d,  uil,  num_points) ;  %  range  of  alpha  in 
degrees 

%  Vavg3  =  ( (Vm*p) / (2*pi) ) * (1  +  cosd(alpha_d  +  lil)); 

%  Vrms3  =  (Vm/sqrt (2 )) *sqrt (( (2+p) /4 )-( (alpha_d*p) /360 )+.. . 

O, 

o 

( (E*p*sind (2*alpha_d) ) / (4*pi) )+( (F*p*cosd (2*alpha_d) ) / (4*pi) ) ) ; 

%  Vrms  ripple3  =  sqrt (Vrms3 . A2  -  Vavg3.A2); 

%  Vpp  ripple3  =  Vm*sind (alpha  d  +  lil); 

%  Vpp  by_Vavg  percent3  =  (Vpp  ripple3  ./  Vavg3) *100; 

%  plot (alpha_d,  Vpp_by_Vavg_percent3 ,  'g-.',  'LineWidth',  2) 
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24-pulse 


o,  o, 
o  o 


p  =  24;  %  p  =  pulse  count;  i.e.,  p  =  6  for  6-pulse  controlled  rectifier 
Vm  =  1  ; 

lil  =  90- (180/p);  %  lil  =  lower  integration  limit  for  start  of  a  pulse 
uil  =  90+ (180/p);  %  uil  =  upper  integration  limit  for  start  of  a  pulse 
A  =  cosd(lil);  %  constant  to  make  math  cleaner/neater 
E  =  cosd (2*lil) ; 

F  =  sind (2*lil) ; 

alpha  bl  d  =  (uil-lil)/2;  %  first  boundary  angle  in  degrees 
alpha  b2  d  =  lil;  %  second  boundary  angle  in  degrees 
num  points  =  1000;  %  number  of  points  for  alpha 

%  code  for  alpha  <  boundary  angle  1 

alpha  d  =  linspace(0,  alpha  bl  d,  num  points);  %  range  of  alpha  in 
degrees 

Vavgl  =  (A*Vm*p*cosd (alpha  d) ) /pi; 

Vrmsl  =  (Vm/sqrt (2) ) *sqrt (1  +  ( (F*p*cosd (2*alpha_d) ) / (2*pi) ) ) ; 

Vrms  ripplel  =  sqrt (Vrmsl . A2  -  Vavgl. A2); 

Vpp  ripplel  =  Vm* ( 1-sind (alpha  d  +  uil)); 

Vpp  by  Vavg  percentl  =  (Vpp  ripplel  .  /  Vavgl) *100; 
plot (alpha  d,  Vpp  by  Vavg  percentl,  'k' ,  'LineWidth' ,  2) 

%  code  for  boundary  angle  1  <  alpha  <  boundary  angle  2 

alpha  d  =  linspace (alpha  bl  d,  alpha  b2  d,  num  points);  %  range  of 

alpha  in  degrees 

Vavg2  =  (A*Vm*p*cosd (alpha  d))/pi; 

Vrms2  =  (Vm/sqrt (2 )) *sqrt ( 1  +  ( (F*p*cosd (2*alpha_d) ) / (2*pi) ) ) ; 

Vrms  ripple2  =  sqrt (Vrms2 . A2  -  Vavg2.A2); 

Vpp  ripple2  =  Vm* (sind (alpha  d  +  lil)  -  sind(alpha  d  +  uil)); 

Vpp  by  Vavg  percent2  =  (Vpp  ripple2  .  /  Vavg2)*100; 

plot (alpha  d,  Vpp  by_Vavg  percent2,  'k--' ,  'LineWidth',  2) 

%  %  code  for  alpha  >  boundary  angle  2 

%  alpha  d  =  linspace (alpha  b2  d,  uil,  num_points) ;  %  range  of  alpha  in 
degrees 

%  Vavg3  =  ( (Vm*p) / (2*pi) ) * (1  +  cosd(alpha_d  +  lil) ) ; 

%  Vrms3  =  (Vm/sqrt (2 )) *sqrt (( (2+p) /4 )-( (alpha_d*p) /360 )+.. . 

Q. 

O 

( (E*p*sind (2*alpha_d) ) / (4*pi) )+( (F*p*cosd (2*alpha_d) ) / (4*pi) ) ) ; 

%  Vrms  ripple3  =  sqrt (Vrms3 . A2  -  Vavg3.A2); 

%  Vpp  ripple3  =  Vm*sind (alpha^d  +  lil) ; 

%  Vpp  by_Vavg  percent3  =  (Vpp  ripple3  ./  Vavg3) *100; 

%  plot (alpha  d,  Vpp  by  Vavg  percent3,  'k-.',  'LineWidth',  2) 

%%  96-pulse 

p  =  96;  %  p  =  pulse  count;  i.e.,  p  =  6  for  6-pulse  controlled  rectifier 
Vm  =  1  ; 

lil  =  90- (180/p);  %  lil  =  lower  integration  limit  for  start  of  a  pulse 
uil  =  90+ (180/p);  %  uil  =  upper  integration  limit  for  start  of  a  pulse 
A  =  cosd(lil);  %  constant  to  make  math  cleaner/neater 
E  =  cosd (2*lil) ; 
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F  =  sind (2*lil) ; 

alpha  bl  d  =  (uil-lil)/2;  %  first  boundary  angle  in  degrees 
alpha  b2  d  =  lil;  %  second  boundary  angle  in  degrees 
num  points  =  1000;  %  number  of  points  for  alpha 

%  code  for  alpha  <  boundary  angle  1 

alpha  d  =  linspace(0,  alpha  bl  d,  num  points);  %  range  of  alpha  in 
degrees 

Vavgl  =  (A*Vm*p*cosd (alpha  d) ) /pi; 

Vrmsl  =  (Vm/sqrt (2) ) *sqrt (1  +  ( (F*p*cosd (2*alpha_d) ) / (2*pi) ) ) ; 

Vrms  ripplel  =  sqrt (Vrmsl . A2  -  Vavgl. A2); 

Vpp  ripplel  =  Vm* ( 1-sind (alpha  d  +  uil)); 

Vpp  by  Vavg  percentl  =  (Vpp  ripplel  .  /  Vavgl) *100; 
plot (alpha  d,  Vpp  by  Vavg  percentl,  'b' ,  'LineWidth' ,  2) 

%  code  for  boundary  angle  1  <  alpha  <  boundary  angle  2 

alpha  d  =  linspace (alpha  bl  d,  alpha  b2  d,  num  points);  %  range  of 

alpha  in  degrees 

Vavg2  =  (A*Vm*p*cosd (alpha  d) ) /pi; 

Vrms2  =  (Vm/sqrt (2 )) *sqrt ( 1  +  ( (F*p*cosd (2*alpha_d) ) / (2*pi) ) ) ; 

Vrms  ripple2  =  sqrt (Vrms2 . A2  -  Vavg2.A2); 

Vpp  ripple2  =  Vm* (sind (alpha  d  +  lil)  -  sind(alpha  d  +  uil)); 

Vpp  by  Vavg  percent2  =  (Vpp  ripple2  .  /  Vavg2)*100; 

plot (alpha  d,  Vpp  by  Vavg  percent2,  'b--' ,  'LineWidth',  2) 

%  %  code  for  alpha  >  boundary  angle  2 

%  alpha  d  =  linspace  (alpha  b2__d,  uil,  num  points);  %  range  of  alpha  in 
degrees 

%  Vavg3  =  ( (Vm*p) / (2*pi) ) * (1  +  cosd(alpha_d  +  lil)); 

%  Vrms3  =  (Vm/sqrt (2 )) *sqrt (( (2+p) /4 )-( (alpha_d*p) /360 )+.. . 

O, 

o 

( (E*p*sind (2*alpha_d) ) / (4*pi) )+( (F*p*cosd (2*alpha_d) ) / (4*pi) ) ) ; 

%  Vrms  ripple3  =  sqrt (Vrms3 . A2  -  Vavg3.A2); 

%  Vpp  ripple3  =  Vm*sind (alpha  d  +  lil); 

%  Vpp  by_Vavg  percent3  =  (Vpp  ripple3  ./  Vavg3) *100; 

%  plot (alpha  d,  Vpp  by  Vavg  percent3,  'b-.',  'LineWidth',  2) 

%%  Plot  Formatting 

grid  on 

axis (  [ 0  30  0  60 ] ) 

%title ( ' %V_p_p  /_a  v  g  %  vs .  \alpha  for  an  n-Pulse  Controlled  Rectifier 

with  a  Purely  Resistive  Load'  ) 

xlabel ( ' \alpha  (degrees)') 

ylabel ( ' %V_p_p_/_a_v_g  (%)  '  ) 

legend (' 6-pulse,  \alpha  \leq  30A\circ' ,  . .  . 

'12-pulse,  \alpha  \leq  15A\circ' ,  '12-pulse,  15A\circ  \leq 
\alpha  \leq  75A\circ'  ,  . .  . 

'18-pulse,  \alpha  \leq  7.5A\circ',  '18-pulse,  7.5A\circ  \leq 
\alpha  \leq  80A\circ' ,  .  .  . 

'24-pulse,  \alpha  \leq  3.75A\circ',  '24-pulse,  3.75A\circ  \leq 
\alpha  \leq  82.5A\circ' ,  .  .  . 

'96-pulse,  \alpha  \leq  1 . 875A\circ' ,  '96-pulse,  1.875A\circ  \leq 
\alpha  \leq  88 . 125A\circ' , . . . 

'location',  'Northwest') 
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%  full  legend  below 

%  legend (' 6-pulse,  \alpha  \leq  30A\circ' ,  '6-pulse,  30A\circ  \leq 
\alpha  \leq  60A\circ'  ,  '6-pulse,  \alpha  \geq  60A\circ' , . . . 

%  '12-pulse,  \alpha  \leq  15A\circ' ,  '12-pulse,  15A\circ  \leq 

\alpha  \leq  75A\circ' ,  '12-pulse,  \alpha  \geq  75A\circ' , . . . 

%  '18-pulse,  \alpha  \leq  7.5A\circ',  '18-pulse,  7.5A\circ  \leq 

\alpha  \leq  80A\circ' ,  '18-pulse,  \alpha  \geq  80A\circ' ,  . . . 

%  '24-pulse,  \alpha  \leq  3.75A\circ',  '24-pulse,  3.75A\circ  \leq 

\alpha  \leq  82.5A\circ',  '24-pulse,  \alpha  \geq  82 . 5 A\circ' , . . . 

%  '96-pulse,  \alpha  \leq  1 . 875A\circ' ,  '96-pulse,  1.875A\circ 

\leq  \alpha  \leq  88 . 125A\circ' ,  '96-pulse,  \alpha  \geq 
88 . 125A\circ' ,  . . . 

%  'location',  'Northwest') 

hold  off 

G.  VOLTAGE  RIPPLE  CALCULATIONS  AND  PLOT  FOR  A  CONTROLLED 
6-PULSE  RECTIFIER  WITH  A  TWO-POLE  LC  OUTPUT  FILTER 


clear 

clc 

hold  on 

%%  initial  constants 

alpha  u  =  30;  %  upper  limit  on  firing  angle  alpha 
num  points  =  10000;  %  number  of  points  for  alpha 
Vm  =  (100/sqrt (2) ) *sqrt (2)  ; 

L0  =  25e-3;  %  inductance  of  load  inductor 
CO  =  50e-6;  %  capacitance  of  load  capacitor 
R0  =  10;  %  resistance  of  load  resistor 

LL  =  L0- ( . 5*L0 ) ; 

LH  =  L0+ ( . 5*L0 ) ; 

CL  =  CO-  ( .5*C0) ; 

CH  =  C0+ ( ,5*C0) ; 

%%  6p_L0_C0 

p  =  6;  %  p  =  pulse  count;  i.e.,  p  =  6  for  6-pulse  controlled  rectifier 
L  =  L0;  %  inductance  of  load  inductor 
C  =  CO;  %  capacitance  of  load  capacitor 
R  =  R0;  %  resistance  of  load  resistor 

lil  =  90- (180/p);  %  lil  =  lower  integration  limit  for  start  of  a  pulse 
uil  =  90+ (180/p);  %  uil  =  upper  integration  limit  for  start  of  a  pulse 
A  =  cosd(lil);  %  constant  to  make  math  cleaner/neater 
F  =  sind (2*lil) ; 

Z_C  =  1  /  ( j *120*p*pi*C) ; 

Z_L  =  j *120*p*pi*L; 

Z_R_p_C  =  (ZJZ  *  R)/(Z_C  +  R)  ; 
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Ztot  =  Z_L  +  Z_R_p_C; 
mag_Ztot  =  abs(Ztot); 

alpha  d  =  linspace(0,  alpha  u,  num  points);  %  range  of  alpha  in  degrees 
Vavg  =  (A*Vm*p*cosd (alpha_d) ) /pi; 

Vrms  =  (Vm/sqrt (2) ) *sqrt (1  +  ( (F*p*cosd (2*alpha_d) ) / (2*pi) ) ) ; 

Vrms  ripple  =  sqrt (Vrms . A2  -  Vavg.A2); 

Vac  =  Vrms  ripple; 

mag  lac  =  Vac  /  mag  Ztot; 

VCrms  =  mag  lac  *  abs(Z  C) ; 

Vpp  ripple  =  VCrms  *  2  *  sqrt (2); 

Vpp_by_Vavg_percent  =  (Vpp_ripple  .  /  Vavg) *100; 

plot (alpha  d,  Vpp  by  Vavg  percent,  'k' ,  'LineWidth' ,  2) 

%%  6p_LL_C0 

p  =  6;  %  p  =  pulse  count;  i.e.,  p  =  6  for  6-pulse  controlled  rectifier 
L  =  LL;  %  inductance  of  load  inductor 
C  =  CO;  %  capacitance  of  load  capacitor 
R  =  R0;  %  resistance  of  load  resistor 

lil  =  90- (180/p);  %  lil  =  lower  integration  limit  for  start  of  a  pulse 
uil  =  90+ (180/p);  %  uil  =  upper  integration  limit  for  start  of  a  pulse 
A  =  cosd(lil);  %  constant  to  make  math  cleaner/neater 
F  =  sind (2*lil) ; 

Z_C  =  1  /  ( j *120*p*pi*C) ; 

Z_L  =  j  *120*p*pi*L; 

Z_R_p_C  =  (Z_C  *  R)/(Z_C  +  R) ; 

Ztot  =  Z_L  +  Z_R_p_C; 
mag_Ztot  =  abs(Ztot); 

alpha  d  =  linspace(0,  alpha  u,  num  points);  %  range  of  alpha  in  degrees 
Vavg  =  (A*Vm*p*cosd (alpha_d) ) /pi; 

Vrms  =  (Vm/sqrt (2 )) *sqrt ( 1  +  ( (F*p*cosd (2*alpha_d) ) / (2*pi) ) ) ; 

Vrms  ripple  =  sqrt (Vrms. A2  -  Vavg.A2); 

Vac  =  Vrms  ripple; 

mag  lac  =  Vac  /  mag  Ztot; 

VCrms  =  mag  lac  *  abs(Z  C) ; 

Vpp  ripple  =  VCrms  *  2  *  sqrt (2); 

Vpp_by_Vavg_percent  =  (Vpp_ripple  .  /  Vavg) *100; 

plot (alpha  d,  Vpp  by_Vavg  percent,  'r' ,  'LineWidth',  2) 

%%  6p_LH_C0 

p  =  6;  %  p  =  pulse  count;  i.e.,  p  =  6  for  6-pulse  controlled  rectifier 
L  =  LH;  %  inductance  of  load  inductor 
C  =  CO;  %  capacitance  of  load  capacitor 
R  =  R0;  %  resistance  of  load  resistor 

lil  =  90- (180/p);  %  lil  =  lower  integration  limit  for  start  of  a  pulse 
uil  =  90+ (180/p);  %  uil  =  upper  integration  limit  for  start  of  a  pulse 
A  =  cosd(lil);  %  constant  to  make  math  cleaner/neater 
F  =  sind (2*lil) ; 
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range  of  alpha  in  degrees 


Z_C  =  1  /  ( j *120*p*pi*C) ; 

Z_L  =  j  *120*p*pi*L; 

Z_R_p_C  =  (Z_C  *  R)/(Z_C  +  R) ; 

Ztot  =  Z_L  +  Z_R_p_C ; 
mag_Ztot  =  abs(Ztot); 

alpha  d  =  linspace(0,  alpha  u,  num  points);  % 

Vavg  =  (A*Vm*p*cosd (alpha_d) ) /pi; 

Vrms  =  (Vm/sqrt (2) ) *sqrt (1  +  ( (F*p*cosd (2*alpha_d) ) / (2*pi) ) ) ; 

Vrms  ripple  =  sqrt(Vrms.A2  -  Vavg.A2); 

Vac  =  Vrms  ripple; 

mag  lac  =  Vac  /  mag  Ztot; 

VCrms  =  mag  lac  *  abs(Z  C) ; 

Vpp  ripple  =  VCrms  *  2  *  sqrt(2); 

Vpp  by  Vavg  percent  =  (Vpp  ripple  .  /  Vavg) *100; 

plot (alpha  d,  Vpp  by  Vavg  percent,  'b' ,  'LineWidth' ,  2) 

%%  6p_L0_CL 

p  =  6;  %  p  =  pulse  count;  i.e.,  p  =  6  for  6-pulse  controlled  rectifier 
L  =  L0;  %  inductance  of  load  inductor 
C  =  CL;  %  capacitance  of  load  capacitor 
R  =  R0;  %  resistance  of  load  resistor 

lil  =  90- (180/p);  %  lil  =  lower  integration  limit  for  start  of  a  pulse 
uil  =  90+ (180/p);  %  uil  =  upper  integration  limit  for  start  of  a  pulse 
A  =  cosd(lil);  %  constant  to  make  math  cleaner/neater 
F  =  sind (2*lil) ; 

Z_C  =  1  /  ( j *120*p*pi*C) ; 

Z_L  =  j  *120*p*pi*L; 

Z_R_p_C  =  (Z_C  *  R)/(Z_C  +  R) ; 

Ztot  =  Z_L  +  Z_R_p_C ; 
mag_Ztot  =  abs(Ztot); 

alpha  d  =  linspace(0,  alpha  u,  num  points);  %  range  of  alpha  in  degrees 
Vavg  =  (A*Vm*p*cosd (alpha_d) ) /pi; 

Vrms  =  (Vm/sqrt (2 )) *sqrt ( 1  +  ( (F*p*cosd (2*alpha_d) ) / (2*pi) ) ) ; 

Vrms  ripple  =  sqrt(Vrms.A2  -  Vavg.A2); 

Vac  =  Vrms  ripple; 

mag  lac  =  Vac  /  mag  Ztot; 

VCrms  =  mag  lac  *  abs(Z  C) ; 

Vpp  ripple  =  VCrms  *  2  *  sqrt(2); 

Vpp  by  Vavg  percent  =  (Vpp  ripple  ./  Vavg) *100; 

plot (alpha  d,  Vpp  by  Vavg  percent,  'k — ' ,  'LineWidth',  2) 

%%  6p_L0_CH 

p  =  6;  %  p  =  pulse  count;  i.e.,  p  =  6  for  6-pulse  controlled  rectifier 
L  =  L0;  %  inductance  of  load  inductor 
C  =  CH;  %  capacitance  of  load  capacitor 
R  =  R0;  %  resistance  of  load  resistor 
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lil  =  90- (180/p);  %  lil  =  lower  integration  limit  for  start  of  a  pulse 
uil  =  90+ (180/p);  %  uil  =  upper  integration  limit  for  start  of  a  pulse 
A  =  cosd(lil);  %  constant  to  make  math  cleaner/neater 
F  =  sind (2*lil) ; 

Z_C  =  1  /  ( j *120*p*pi*C) ; 

Z_L  =  j *120*p*pi*L; 

Z_R_p_C  =  (Z_C  *  R)/(Z_C  +  R) ; 

Ztot  =  Z_L  +  Z_R_p_C; 
mag_Ztot  =  abs(Ztot); 

alpha  d  =  linspace(0,  alpha  u,  num  points);  %  range  of  alpha  in  degrees 
Vavg  =  (A*Vm*p*cosd (alpha_d) ) /pi; 

Vrms  =  (Vm/sqrt (2) ) *sqrt (1  +  ( (F*p*cosd (2*alpha_d) ) / (2*pi) ) ) ; 

Vrms  ripple  =  sqrt(Vrms.A2  -  Vavg.A2); 

Vac  =  Vrms  ripple; 
mag_Iac  =  Vac  /  mag_Ztot; 

VCrms  =  mag  lac  *  abs(Z  C) ; 

Vpp  ripple  =  VCrms  *  2  *  sqrt(2); 

Vpp  by  Vavg  percent  =  (Vpp  ripple  ./  Vavg) *100; 

plot (alpha  d,  Vpp  by  Vavg  percent,  'm'  ,  'LineWidth' ,  2) 

%%  Plot  Formatting 


grid  on 

axis (  [0  30  0  20] ) 

%title ( 'Vpp-ripple/Vavg  %  vs.  \alpha  for  a  6-Pulse  Controlled  Rectifier 
with  a  Two-Pole  LC  Output  Filter'  ) 
xlabel ( ' \alpha  (degrees)') 
ylabel ( '%V_p_p_/_V_a_v_g  (%) ' ) 


legend ( 'L 

0 

c 

0' ,  .  .  . 

'L 

L 

c 

0' ,  .  .  . 

'L 

H 

c 

0' ,  .  .  . 

'L 

0 

c 

L' ,  .  .  . 

'L 

0 

c 

H' ,  ■  ■  ■ 

'location',  'Northwest') 
hold  off 

H.  VOLTAGE  RIPPLE  CALCULATIONS  AND  PLOT  FOR  A  CONTROLLED 
6-,  12- ,  18- ,  24-,  AND  96-PULSE  RECTIFIER  WITH  A  TWO- 

POLE  LC  OUTPUT  FILTER 


clear 

clc 

hold  on 

%%  initial  constants 

alpha  u  =  30;  %  upper  limit  on  firing  angle  alpha 
num  points  =  10000;  %  number  of  points  for  alpha 
Vm  =  (100/sqrt (2) ) *sqrt (2) ; 
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LO  =  25e-3;  %  inductance  of  load  inductor 
CO  =  50e-6;  %  capacitance  of  load  capacitor 
RO  =  10;  %  resistance  of  load  resistor 

%%  6p_L0_C0 

p  =  6;  %  p  =  pulse  count;  i.e.,  p  =  6  for  6-pulse  controlled  rectifier 
L  =  LO;  %  inductance  of  load  inductor 
C  =  CO;  %  capacitance  of  load  capacitor 
R  =  RO;  %  resistance  of  load  resistor 

lil  =  90-(180/p);  %  lil  =  lower  integration  limit  for  start  of  a  pulse 
uil  =  90+ (180/p);  %  uil  =  upper  integration  limit  for  start  of  a  pulse 
A  =  cosd(lil);  %  constant  to  make  math  cleaner/neater 
F  =  sind (2*lil)  ; 

Z_C  =  1  /  ( j *120*p*pi*C) ; 

Z_L  =  j  *120*p*pi*L; 

Z_R_p_C  =  (Z_C  *  R)/(Z_C  +  R) ; 

Ztot  =  Z_L  +  Z_R_p_C; 
mag_Ztot  =  abs(Ztot); 

alpha  d  =  linspace(0,  alpha  u,  num  points);  %  range  of  alpha  in  degrees 
Vavg  =  (A*Vm*p*cosd (alpha_d) ) /pi; 

Vrms  =  (Vm/sqrt (2) ) *sqrt (1  +  ( (F*p*cosd (2*alpha_d) ) / (2*pi) ) ) ; 

Vrms  ripple  =  sqrt (Vrms . A2  -  Vavg.A2); 

Vac  =  Vrms  ripple; 
mag_Iac  =  Vac  /  mag_Ztot; 

VCrms  =  mag  lac  *  abs(Z  C) ; 

Vpp  ripple  =  VCrms  *  2  *  sqrt (2); 

Vpp  by  Vavg  percent  =  (Vpp  ripple  .  /  Vavg) *100; 

plot (alpha  d,  Vpp  by  Vavg  percent,  'r' ,  'LineWidth' ,  2) 

%%  12pJL0_C0 

p  =  12;  %  p  =  pulse  count;  i.e.,  p  =  6  for  6-pulse  controlled  rectifier 
L  =  LO;  %  inductance  of  load  inductor 
C  =  CO;  %  capacitance  of  load  capacitor 
R  =  RO;  %  resistance  of  load  resistor 

lil  =  90-(180/p);  %  lil  =  lower  integration  limit  for  start  of  a  pulse 
uil  =  90+(180/p);  %  uil  =  upper  integration  limit  for  start  of  a  pulse 
A  =  cosd(lil);  %  constant  to  make  math  cleaner/neater 
F  =  sind (2*lil) ; 

Z_C  =  1  /  ( j *120*p*pi*C) ; 

Z_L  =  j  *120*p*pi*L; 

Z_R_p_C  =  (ZJZ  *  R)/(Z_C  +  R)  ; 

Ztot  =  Z_L  +  Z_R_p_C; 
mag_Ztot  =  abs(Ztot); 

alpha  d  =  linspace(0,  alpha  u,  num  points);  %  range  of  alpha  in  degrees 
Vavg  =  (A*Vm*p*cosd (alpha_d) ) /pi; 

Vrms  =  (Vm/sqrt (2 )) *sqrt ( 1  +  ( (F*p*cosd (2*alpha_d) ) / (2*pi) ) ) ; 
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Vrms  ripple  =  sqrt(Vrms.A2  -  Vavg.A2); 

Vac  =  Vrms  ripple; 

mag  lac  =  Vac  /  mag  Ztot; 

VCrms  =  mag  lac  *  abs(Z  C) ; 

Vpp  ripple  =  VCrms  *  2  *  sqrt(2); 

Vpp_by_Vavg_percent  =  (Vpp_ripple  .  /  Vavg) *100; 

plot (alpha  d,  Vpp  by_Vavg  percent,  'm'  ,  'LineWidth' ,  2) 

%%  18pJL0_C0 

p  =  18;  %  p  =  pulse  count;  i.e.,  p  =  6  for  6-pulse  controlled  rectifier 
L  =  LO;  %  inductance  of  load  inductor 
C  =  CO;  %  capacitance  of  load  capacitor 
R  =  RO;  %  resistance  of  load  resistor 

lil  =  90- (180/p);  %  lil  =  lower  integration  limit  for  start  of  a  pulse 
uil  =  90+ (180/p);  %  uil  =  upper  integration  limit  for  start  of  a  pulse 
A  =  cosd(lil);  %  constant  to  make  math  cleaner/neater 
F  =  sind (2*lil)  ; 

Z_C  =  1  /  ( j *120*p*pi*C) ; 

Z_L  =  j *120*p*pi*L; 

Z_R_p_C  =  (ZJZ  *  R)/(Z_C  +  R)  ; 

Ztot  =  Z_L  +  Z_R_p_C; 
mag_Ztot  =  abs(Ztot); 

alpha  d  =  linspace(0,  alpha  u,  num  points);  %  range  of  alpha  in  degrees 
Vavg  =  (A*Vm*p*cosd (alpha_d) ) /pi; 

Vrms  =  (Vm/sqrt (2) ) *sqrt (1  +  ( (F*p*cosd (2*alpha_d) ) / (2*pi) ) ) ; 

Vrms  ripple  =  sqrt(Vrms.A2  -  Vavg.A2); 

Vac  =  Vrms  ripple; 

mag  lac  =  Vac  /  mag  Ztot; 

VCrms  =  mag  lac  *  abs(Z  C) ; 

Vpp  ripple  =  VCrms  *  2  *  sqrt(2); 

Vpp_by_Vavg_percent  =  (Vpp_ripple  ./  Vavg) *100; 

plot (alpha  d,  Vpp  by_Vavg  percent,  'g' ,  'LineWidth',  2) 

%%  24p_L0_C0 

p  =  24;  %  p  =  pulse  count;  i.e.,  p  =  6  for  6-pulse  controlled  rectifier 
L  =  LO;  %  inductance  of  load  inductor 
C  =  CO;  %  capacitance  of  load  capacitor 
R  =  RO;  %  resistance  of  load  resistor 

lil  =  90— (180/p) ;  %  lil  =  lower  integration  limit  for  start  of  a  pulse 
uil  =  90+ (180/p);  %  uil  =  upper  integration  limit  for  start  of  a  pulse 
A  =  cosd(lil);  %  constant  to  make  math  cleaner/neater 
F  =  sind (2*lil)  ; 

Z_C  =  1  /  ( j*120*p*pi*C) ; 

Z_L  =  j *120*p*pi*L; 

Z_R_p_C  =  (Z_C  *  R)/(Z_C  +  R) ; 

Ztot  =  Z_L  +  Z_R_p_C; 
mag_Ztot  =  abs(Ztot); 
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alpha  d  =  linspace(0,  alpha  u,  num  points);  %  range  of  alpha  in  degrees 
Vavg  =  (A*Vm*p*cosd (alpha_d) ) /pi; 

Vrms  =  (Vm/sqrt (2) ) *sqrt (1  +  ( (F*p*cosd (2*alpha_d) ) / (2*pi) ) ) ; 

Vrms  ripple  =  sqrt(Vrms.A2  -  Vavg.A2); 

Vac  =  Vrms  ripple; 

mag  lac  =  Vac  /  mag  Ztot; 

VCrms  =  mag  lac  *  abs(Z  C) ; 

Vpp  ripple  =  VCrms  *  2  *  sqrt(2); 

Vpp  by  Vavg  percent  =  (Vpp  ripple  ./  Vavg) *100; 

plot (alpha  d,  Vpp  by_Vavg  percent,  'k' ,  'LineWidth' ,  2) 

%%  96p_L0_C0 

p  =  96;  %  p  =  pulse  count;  i.e.,  p  =  6  for  6-pulse  controlled  rectifier 
L  =  L0;  %  inductance  of  load  inductor 
C  =  CO;  %  capacitance  of  load  capacitor 
R  =  R0;  %  resistance  of  load  resistor 

lil  =  90- (180/p);  %  lil  =  lower  integration  limit  for  start  of  a  pulse 
uil  =  90+ (180/p);  %  uil  =  upper  integration  limit  for  start  of  a  pulse 
A  =  cosd(lil);  %  constant  to  make  math  cleaner/neater 
F  =  sind (2*lil)  ; 

Z_C  =  1  /  ( j*120*p*pi*C) ; 

Z_L  =  j *120*p*pi*L; 

Z_R_p_C  =  (Z_C  *  R)/(Z_C  +  R) ; 

Ztot  =  Z_L  +  Z_R_p_C ; 
mag_Ztot  =  abs(Ztot); 

alpha  d  =  linspace(0,  alpha  u,  num  points);  %  range  of  alpha  in  degrees 
Vavg  =  (A*Vm*p*cosd (alpha_d) ) /pi; 

Vrms  =  (Vm/sqrt (2 )) *sqrt ( 1  +  ( (F*p*cosd (2*alpha_d) ) / (2*pi) ) ) ; 

Vrms  ripple  =  sqrt(Vrms.A2  -  Vavg.A2); 

Vac  =  Vrms  ripple; 

mag  lac  =  Vac  /  mag  Ztot; 

VCrms  =  mag  lac  *  abs(Z  C) ; 

Vpp  ripple  =  VCrms  *  2  *  sqrt(2); 

Vpp  by  Vavg  percent  =  (Vpp  ripple  ./  Vavg) *100; 

plot (alpha  d,  Vpp  by  Vavg  percent,  'b' ,  'LineWidth',  2) 

%%  Plot  Formatting 


grid  on 

axis (  [0  30  0  9] ) 

%title ( 'Vpp-ripple/Vavg  %  vs.  \alpha  for  a  6-Pulse  Controlled  Rectifier 

with  a  Two-Pole  LC  Output  Filter' ) 

xlabel (  ' \alpha  (degrees)') 

ylabel (  ' %V_p_p_/_V_a_v_g  (%)  '  ) 

legend ( '6-pulse'  ,  .  .  . 

'12-pulse' ,  .  .  . 

'18-pulse' ,  .  .  . 

'24-pulse' ,  .  .  . 
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'96-pulse' , 
'location' , 
hold  off 


'Northwest' ) 
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